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ABSTRACT 

We discuss three new equations of state (EOS) in core-collapse supernova simulations. The new EOS 
are based on the nuclear statistical equilibrium model of Hempel and Schaffner-Bielich (HS), which 
includes excluded volume effects and relativistic mean-field (RMF) interactions. We consider the RMF 
parameterizations TM1, TMA, and FSUgold. These EOS are implemented into our spherically sym- 
metric core-collapse supernova model, which is based on general relativistic radiation hydrodynamics 
and three-flavor Boltzmann neutrino transport. The results obtained for the new EOS are compared 
with the widely used EOS of H. Shen et al. and Lattimer & Swesty. The systematic comparison shows 
that the model description of inhomogeneous nuclear matter is as important as the parameterization 
of the nuclear interactions for the supernova dynamics and the neutrino signal. Furthermore, several 
new aspects of nuclear physics are investigated: the HS EOS contains distributions of nuclei, includ- 
ing nuclear shell effects. The appearance of light nuclei, e.g., deuterium and tritium is also explored, 
which can become as abundant as alphas and free protons. In addition, we investigate the black hole 
formation in failed core-collapse supernovae, which is mainly determined by the high-density EOS. We 
find that temperature effects lead to a systematically faster collapse for the non-relativistic LS EOS 
in comparison to the RMF EOS. We deduce a new correlation for the time until black hole formation, 
which allows to determine the maximum mass of proto-neutron stars, if the neutrino signal from such 
a failed supernova would be measured in the future. This would give a constraint for the nuclear EOS 
at finite entropy, complementary to observations of cold neutron stars. 



1. INTRODUCTION 

Supernova explosions of stars more massive than 8 M Q 
are an active subject of research in astrophysics. The 
core-collapse supernova problem is related to the revival 
of the stalled bounce shock. It forms when the collapsing 
stellar core bounces back above normal nuclear matter 
density. Explosions in spherical symmetry have been ob- 
tained on ly for the low mass 8. 8 M Pl O-Ne-Mg core from 
Nomotol (U98l H98l ll~987ft by Kitaura eTall pOOl) and 
Fischer et al.l ([2010). More massive progenitor stars have 
extended high-density Si-layers surrounding the central 
iron core. The post-bounce evolution leads to a mass ac- 
cretion period that lasts for several hundred milliseconds. 
For such iron-core progenitor stars explosions do not oc- 
cur in spherically symmetric simulations. The bounce 
shock continuously looses energy by dissociation of heavy 
nuclei and neutrino heating is not sufficient enough to re- 
vive the standing bounce shock. 

In addition to the s tandard scenario of n eutrino- 
driven explosions (see iBethe &: Wilson! (|1985ft ). sev- 
eral alternative explosion mechanisms have been pro- 
po sed. These are the magne to-rotational mechanism 
by iLeBlanc and Wilson! (119701) and the acoustic mech- 
anism by IBurrows et al.l (|2006ft . All of which are work- 
ing in multiple spatial dimensions. Multi-dimensional 
core-collapse supernova models, based on sophisticated 
neutrino transport approximations and general relativ- 
ity, have become available only recently. These models 
have shown to increase t he ne utr ino heating e f ficienc y 
(see, e.g., iMiller et all (11993ft. iHerant et all (|1994ft . 
Uanka fc Muellerl (|1996l )" IBurrows et all (|2007l )) and help 



to aid the understanding of asp h erical explosions ( see, fo r 
exa mple. iBruenn et al (12006ft: iMarek fc Jankal (|2009ft : 
and lMiiller et al.l (|2010ft ). 

Apart from the dimensionality, the simulations are also 
affected by uncertainties in the nuclear physics involved, 
which can be divided in weak processes and the equa- 
tion of state (EOS). Changes in the nuclear physics input 
can have dramatic consequences and can even generate 
explosions in sphe rically symm e tric m odel s as was illus- 
trated recently bv lSagert et al.l (2009) and lFischer et all 
(|2011l) . The authors explored the state of matter around 
and above nuclear matter density and investigated the 
possibility of the quark-hadron phase transition during 
the early post-bounce evolution of low and intermediate 
mass iron-core progenitor stars. It was found that the 
phase transition leads to the formation of a strong hy- 
drodynamic shock wave which triggers the explosion. In 
the present article we study new aspects of the hadronic 
EOS and their consequences for the supernova dynamics. 

Until today, there are only a few hadronic EOS avail- 
able which cover a sufficiently large domain in density, 
temperature, and electron fraction so that they can be 
used in simulations of core-collapse supernovae. Usu- 
ally, in such simulations the EOS is implemented in 
form of a table. The most prominent among th e m are 
the two "clas sic" EOS from iLattimer fc Swestvl (|1991ft 
(LS) and H. IShen et all d!998a|)bj) (STOS) which are 
commonly used in computational astrophysics. The LS 
EOS is provided in form of programming routines for 
three different values of the nuclear incompressibility K 
of 180, 220, and 375 MeV. STOS is an EOS in tabu- 
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lar form which was calc ulated for the relativisti c mean- 
field (RMF) model TM1 (jSugahara fc TokilU99l . There 
has been a lot of progress for the supernova EOS in 
the l ast years. Recently , the new hadronic EOS tables 
of G. iShen et all (|2011bl |ah became available, which are 
based on the virial expansion and two different RMF in- 
teractions in the Hartree app roximation. In the ir first 
article they applied NL3 (G. IShen et al.l (l2011bl) N >. fo l- 
lowed by a table for FSUeold (G. IShen et al.l (1201 lal) ). 
Furthermore also an update of the original STOS EOS 
table appeared which has improved res olution, accuracy 
and grid-spacing (H. IShen et aD fj201 lh ) . In this article 
also another table was presented in which lambda hy- 
perons are taken into account. Similarly, several other 
extensions of the STOS tab le are available, eith er with 
the formation of hyperons (Ishizuka et al.ll2008fl or the 
inclus ion of the QCD ph ase transition fNakazat o et al.1 
[20081: ISagert et all [20091: iFischer et all miTW . In the 
present investigation we apply th e new EOS tables of 
IHempel fe Schaffner-Biehchl ((2010l) . 

An important aspect of the supernova EOS is the 
formation of light nuclei and their properties in the 
hot and dense medium. Note that the two clas- 
sic supernova EOS include only alpha particles of all 
possible light nuclei, w hich are imp l ement ed with ex- 
cluded volume effects. iTypel et al.l (|2010T ) use a gen- 
eraliz ed RMF and a quantum statistical model {Ropkc 
2011) to study supernova matter taking into account 
the most important light nuclei with mass number 
A < 4. Significant differences in the composition 
were foun d compared to STOS. simila r as in the recent 
studie s of iHorowitz fc Schwenkl (120061) . lO'Connor et al.l 
(120071). iSumivoshi fc Ropkd (120081) . iHeckel et al.l (120091). 
and IHempel fc Schaffner-Bielichl (|2010t ). Also thermo- 
dynamic variables, like ,e.g., the symmetry energy, are 
modified due to the appearance of light nuclei. A com- 
parison between the excluded volume approach with all 
light nuclei, which is also used in th e HS E O S, an d the 
quantum m a ny-bo dy approac hes of iRopkd (120111) and 
ITypel et al.l (|2010f ) is given bv IHempel et all (|2011h . It 
is found that the dissolution of light nuclei happens at 
similar densities and that the excluded volume approach 
gives an acceptable representation of the real microscopic 
medium effects, like, e.g., Pauli-blocking. Recently, the 
importance of light nuclei in supernova matter was also 
shown by a he avy-ion collision e xperi ment. In the isoscal- 
ing analysis of lNatowitz et all (|2010h the measured sym- 
metry energy can only be explained, if the formation of 
light nuclei is taken properly into account. Quite re- 
markably, the deduced temperatures and densities cor- 
respond to conditions which are typical for matter in 
core-collapse supernovae. In the supernova environment 
light nuclei can possibly influence the neutrino transport 
and consequently the supernova neutrin o signal and dy- 
namic s , see, e.g.. |0'Connor et al. (|2007). lArcones et al 
(2008), S umivoshi fc Ropkd (12008 s ).' and lNakamura et al 
2009). However, so far there are no investigations of 



core-collapse supernovae, which consistently take into ac- 
count all light nuclei. With this article we further ad- 
vance in this directio n, by applying the EOS table s of the 
statistical model of IHempel "fc S chaffner -Bielichl ([2010D 
(HS) to simulations of massive stars. 

Statistical models are characterized by having a 
distribution of nuclei, in which light nuclei can also be 



included. On the other hand, the complexity and variety 
of the involved aspects of nuclear physics require some 
modeling. Important aspects in the supernova context 
are for example the interactions of unbound nucleons, 
excluded volume effects, excited states, Coulomb inter- 
actions, nuclear binding energies, shell effects or surfac e 
modifications of nucl ei, see, e. g., Ilshizuka et al.l ([ 2003). 
IBotvina fc Mishustinl (120041) . INadvozhin fc Yudinl 
(12004D. INadvozhin fc Yudinl (120051) . ISouza et al.l (12 009) , 
Blinnikov et al.l (1201 ID. IHempel fc Schaffner-B ielich 
(|2010[) IBotvina fc Mishustinl (I2010D. lArcones et al.l 
(l27ilfl).lRaduta fc Gulminellil (|2010f ). lYudinl (I201l l). and 
iFurusawa et al.l (|2011[ ). Interestingly, similar statistical 
models are also u sed for the analysis of multifragmenta- 
tion experiments dKoonin fc Randrup|[l987t IGrosslll990l: 
IBondorf et aTlll995D. In this article w e appl y the new 
EOS from IHempel fc Schaffner-Bielichl ((2010). The HS 
EOS uses nuclear statistical equilibrium (NSE) with 
excluded volume effects for distributions of nuclei and 
RMF interactions of unbound nucleons. In this model 
all possible light nuclei (e.g., deuterium and tritium) 
and heavy nuclei up to mass number A ~ 330 from 
the neutron to the proton dripline are included. The 
knowledge of the detailed composition makes it possible 
to analyze the impact of different baryon contributions 
during the different phases of core-collapse supernovae. 

The nuclear distributions are especially relevant for 
weak reactions. For example the new treatments o f elec- 
tro n captures on hea vy nuclei bv lLanganke et ahl (|2003l ) 
and lHix et al.l (|2003|) and inelastic neutrin o-nucleon (nu- 
clei) scattering bv lLanganke et~aLl (|2008l ) are all based 
on NSE distributions of nuclei. Idealistically, the nuclear 
physics input in the weak reactions should be as con- 
sistent as possible with the EOS and therefore the two 
parts should be built from the same degrees of freedom. 
However, also in the present study we use additional sim- 
plifications, which will be specified below. 

To classify the characteristic features of HS, we take 
the two aforementioned classic supernova EOS, i.e., LS 
and STOS, as standard references. These two EOS both 
use the single nucleus approximation (SNA), i.e., the dis- 
tribution of he avy nuclei is replaced by a sin gle represen- 
tative nucleus. iBurrows fc Lattimerl fl984D showed that 
the SNA has only a minor effect on thermodynamic vari- 
ables. However, using the SNA instead of nuclear distri- 
butions represents a significant differe nce for the weak re - 
actions, as mentioned already above. ISouza et al.l ([2009) 
found that the SNA gives also a systematically larger 
representative heavy nucleus compared to the average 
of the distribution. There are further limitations of the 
description of nuclei in the two EOS models: STOS is 
based on the Thomas-Fermi and local-density approxi- 
mation and a minimization of the free energy for param- 
eterized nucleon-density profiles within an RMF model. 
LS is based on a non-relativistic liquid drop description 
including surface effects. Both models do not include 
nuclear shell effects and only give an approximate de- 
scription of the iron-group nuclei, which appear at low 
temperatures. This leads to artificial shifts in the EOS 
at the transition to non-NSE. Contrary, due to the use of 
experimental nuclear masses, the HS model reproduces 
the thermodynamic state of the ideal gas of, e.g., 56 Fe 
or 56 Ni (depending on the proton-to-baryon ratio) at 
low temperatures by construction and naturally includes 
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shell effects. Thus, it represents a major improvement 
to other commonly used EOS as it gives a more detailed 
description of the chemical composition and the nuclear 
effects involved at low densities below ~ fO 12 g/cm 3 . 

However, with statistical models it is difficult to de- 
scribe the transition to uniform nuclear matter which 
occurs around half times saturation density, i.e., around 
10 14 g/cm 3 . Microscopic SNA models allow a more de- 
tailed description of the effects which occur at such large 
densities, like the formation of (non-spherical) inhomoge- 
neo us structures, the so-ca lled nuclear pasta phases, see , 
e. g..lRavenhall et alJ(l!983f> andlNewton fc Stone] (|2009f ). 
In lHempel fe Schaffner-Bielichl (201 0) it was shown that 
the thermodynamic variables in the transition region cal- 
culated with the HS model are in satisfactory agreement 
with the results of STOS and LS. From this point of 
view, the present investigation serves as a practical test 
whether the HS EOS is suitable to be used for all possible 
conditions in core-collapse supernova simulations. 

Recently, five new supernova EOS tables were calcu- 
lated with the HS model for different nuclear interac- 
tions and are available onlineQ One of the new HS 
EOS tables is based on the same parameterization TMf 
(jSugahara fc Tokil ["l994) of RMF interactions of the nu- 
cleons as STOS. This makes it possible to investigate 
the impact of the model assumptions for the descrip- 
tion of inhomogeneous nuclear matter below saturation 
density, which is one of the main topics of the present 
study. In addition to TMf and the general features of 
the HS model, we also investigate two othe r RMF pa- 
ramet erizations of nuc lear interactions, TMA dToki et al.l 
(fl995l)) and FSUgo ld (|Todd-Rutel fc Piekarewiczl(j2005l) ; 
iPiekarewiczl (|2007f )). We note that these two interac- 
tions have not been applied in supernova-simulations so 
far. The different nuclear interactions become significant 
at high densities and we analyze to which extent they 
can influence the supernova dynamics and the observ- 
able neutrino signal. For the distinction of the different 
HS EOS we append the used parameterization in brack- 
ets, e.g., "HS (TMA)". To be explicit, sometimes we also 
use "STOS (TM1)" although the STOS table exists only 
for TM1. For the LS EOS we append the used incom- 
pressibility in brackets, e.g., "LS (180)." 

The manuscript is organized as follows. In |J2J we 
present our core-collapse model, including the new HS 
EOS. The results from simulations of the collapse, 
bounce, and the early post-bounce phases of a 15 M© 
progenitor are presented in [J3] We compare HS (TM1) 
with the LS (180) and the STOS (TM1) EOS and also an- 
alyze the neutrino signal. In JQJ we discuss the neutrino 
signal of a 40 M Q progenitor. For this massive progenitor 
we apply five different EOS: LS (180), LS (220), STOS 
(TM1), HS (FSUgold), and HS (TMA). We focus on the 
time until black hole formation, and investigate its rela- 
tion to characteristic properties of the different nuclear 
interactions, before we close with a summary in ij5] 

2. THEORETICAL AND NUMERICAL SETUP 

2.1. Core- collapse Supernova Model 

Our core-collapse model was originally constructed 
based on Newtonian radiation hydrodynamics and three 

1 See http : //phys -merger .physik .unibas . ch/~hempel/eos . html. 



TABLE 1 

Neutrino Reactions Considered, Including References 



Weak Process a 



References 



u e + n <= p + e 
Pe + p ^ n + e + 
v e + (A) <=* (A) + e~ 
v + N/(A)/a f± u' + N/(A) /a 




u + D + N + N^N + N 



Bruenn (19851 : 
Mczzacappa & Bruenn (1993a) 

Bruenn (19851 ; 
Mczzacappa & Bruenn (1993c) 

Bruenn (19851 ; 
Mczzacappa & Mcsscr (1999) 

O. E. B. MesserJ^ 

Buras et al. (20031 ; 
Fischer et al. (20091 



Notes: 

a u = {^ e ,P e ,^/ T ,P M / T } and N = {n,p}. 

b Privatc communications, based on Hanncstad & Raffclt (1998). 

flavor Boltzmann neutrino t ransport (for details see 
iMezzacappa fc Bruennl 1993a, b.c]). It was extended to 
solve the general relativistic equatio ns for both hydro- 
dynam ics and neutrino transport in ILiebendorfer et ID 
(|2001| ). based on the following line-element 



ds 2 = a 2 dt 2 



-] da 2 



r 2 dfl 



in co-moving coordinates system time t and enclosed 
baryon mass a, with metric coefficients lapse function 
a(t,a) and T(t, a), where r' — dr/da. r 2 dfl describes a 
2-sphcrc, in spherical coordinates (9,4>), of radius r(t,a). 
Equations for momentum and energy conservation are 
obtained via the co- variant derivative of the stress-energy 
tensor, V, T ik = 0, with 

T u = p(l + e + J), T aa = p + pK, T at = T ta = pH, 
T 6 ° = T**= p+ 'ip{j-K), 

where p and p are rest-mass density (further discussed 
below) and matter pressure, e is the specific internal 
matter energy with respect to the rest-mass density. J 
and K are the neutrino contributions to the specific in- 
ternal energy and specific pressure. H is the specific 
energy exchange between matter and the neutrino radia- 
tion field. H, J, and K are the neutrino moments and are 
given by momentum integrals of the neutrino distribu - 
tion functions (for details see ILiebendorfer et alj .2004). 
Furthermore, special attention has been devoted to accu- 
ra tely conserve energy, mom entum, and lepton number 
in ILiebendorfer ct al. 



Neutrino contributions, more precisely the evolution 
of the neutrino distribution functions /„, are computed 
solving the Boltzmann transport equation for ultrarel- 
ativistic and mass-less fermions, for each flavor v — 
v e ,i? e ,v ll / T ,i' l j / i T . In addition to the space-time coordi- 
nates, it depends on the momentum angle 9 and energy 
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E. The left-hand side of the Boltzmann equation deter- 
mines the phase-space derivative due to transport, and 
the right-hand side takes collisions into account. The lat- 
ter relates to the weak processes considered, which are 
listed in Table [T] Note that for the alpha particles only 
neutral-current elastic scattering reactions are included. 
The emission of (/x/r)-ncutrino pairs via the annihila- 
tion of t rapped (v e , u e ) has b een implemented into our 
model in Fisc her et al.l (|2009f) . In the Boltzmann trans- 
port representation, for scattering processes we employ 
reaction kernels, e.g., they take into account final state 
blocking and depend on incoming and outgoing neutrino 
angle and energy. For charged-current processes, we em- 
ploy reaction rates which depend only on the incoming 
neutrino energy. 

The radiation hydrodynamics equations, which arc 
given explicitly in the co-moving frame of reference, are 
discretized numerically based on a Lagrangian descrip- 
tion and an adaptive mesh. They are solved implic- 
itly, where the advection scheme used is second-order 
total variation diminishin g. Details can be obtained i n 
ILiebendorfer et al.l (|2002fl and ILiebendorfer et al.l (|2004D . 
For the current study, we include 104 radial adaptive 
mass zones. Resolution is adapted according to high 
gradients in density, velocity, and entropy. Further- 
more, the outer boundary condition is set by applying 
a Schwarzschild metric at the outermost mass zone, as- 
suming vacuum outside. In addition, neutrino transport 
is discretized in 6 angular bins (cos 9 £ { — 1, +1}) and 
20 e nergy bins (E = 3, . . , 300 MeV). 

In ILiebendorf er et al.l (|2004D . the authors compare the 
model with S. Bruenn's multi-group flux limited diffusion 
neutrino transport approximation. A similar compari- 
son with the Boltzmann tran sport used by the Garchin g 
group has been published in ILiebendorfer et al.l {2005). 
These two studies show a good qualitative agreement of 
the spherically symmetric models used by the different 
groups. 

2.2. Use of the EOS 

The EOS in supernova simulations is part of the nu- 
clear physics input, next to weak interactions in the 
neutrino transport. There are two intrinsically differ- 
ent regimes. In NSE, the destruction and production of 
nuclei is in thermal and chemical equilibrium regarding 
strong and electromagnetic interactions. The conditions 
for NSE are achieved typically above a temperature of 
T ~ 0.5 MeV0 For temperatures below T = 0.44 MeV, 
where NSE cannot be applied, we use an ideal Si-gas. 
For the simulations with the HS and STOS EOS in the 
non-NSE regime we calculate the baryon EOS as an ideal 
Maxwell-Boltzmann gas of silicon with the measured nu- 
clear binding energy from lAudi et~aTI pOOl . For LS this 
is not necessary, because the LS routines used in the su- 
pernova model also provide the EOS of the silicon gas. 
The LS routines also handle intrinsically the transition to 
the NSE regime at temperatures of T > 0.44 MeV. For 
the HS and STOS EOS the baryon component in NSE 
is given in tabular form. It is a good approximation to 
treat electrons (and positrons) as a rigid uniform back- 
ground. Thus they can be added as non-interacting ideal 

2 Throughout the article we use natural units, i.e., h = c =ke = 
1, for quantities where it is appropriate. 



Fermi-Dirac gases based on iTimmes fc Arnettl ([1999). 
The blackbody contribution of photons is also taken into 
account. 

The baryon EOS tables are stored for given temper- 
ature T, baryon number density Ub, and total proton 
fraction Y p = n* ot /ns. Because of charge neutrality, in 
the absence of muons the number density of electrons n e 
is given by n e = Y p ub which connects the baryon and 
the electron EOS. Accordingly, Y p is equivalent to the 
electron fraction Y e — n e /riB — Y p . All particle densities 
and fractions used in the article correspond to net densi- 
ties, i.e., the difference between particle and antiparticle 
densities, e.g., n e — n e - — n e +, Y e = Y e - — Y e +. In astro- 
physics, very often instead of the baryon number density 
ris a baryon mass density p is used. This is completely 
equivalent, if p is defined as p = n^m, with an arbitrary 
but constant mass m. We will also use such a "baryon 
mass density" with the usual convention p — riBm u , with 
the atomic mass unit m u = 931.49432 MeV. However, 
we want to stress that p defined in such a way is not 
the rest-mass density, because, e.g., the dominant parti- 
cle species in NSE depends on the thermodynamic state 
of the system. In a relativistic description of the EOS 
with interactions there is no conservation of rest mass, 
but only of the total energy and baryon number. Anal- 
ogously to the baryon mass density one also defines a 
baryonic mass M = Nbttl u which is just a redefinition of 
the total baryon number Nb of the investigated system. 

The HS EOS contains the full distribution of all avail- 
able nuclei, as will be discussed in more detail in the next 
subsection. However, regarding weak processes we only 
consider an average heavy nucleus and an average light 
nucleus. We separate the distribution of all nuclei into 
light nuclei and heavy nuclei by the charge number six, 
i.e. carbon: 

X a = ^ An A , z /n Bl (1) 

A>2,Z<5 

X A = An A ,zln B , (2) 

A>2,Z>6 

with ua,z denoting the number density of a nucleus with 
mass number A and charge number so that the sum 
of the mass fractions of neutrons X n , protons X p , light 
nuclei X a and heavy nuclei X A is unity. The average 
mass and charge of the light and heavy nuclei are: 

( a ) = Y An AtZ / Y n AtZ , (3) 

A>2,Z<5 A>2,Z<5 

( z ) = Z n A,Z / X! nA > z -> ( 4 ) 

A>2,Z<5 A>2,Z<5 

(A) = A n A>z I Y nA ' Z > ( 5 ) 

A>2,Z>6 A>2,Z>6 
( Z ) = Y Z UA - Z 

I ]T n A ,z. (6) 

A>2,Z>6 A>2,Z>6 

We further simplify the neutrino reactions with light 
nuclei by treating all of them as alpha particles. Note 
that for the alpha particles, and thus for all light nu- 
clei, in the present investigation only neutral-current re- 
actions are considered, as listed in Table [TJ Charged- 
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current and other break-up reactions of the weakly bound 
light nuclei may a lso be importan t in su pernova sim- 
ulations, see, e.g., iNakamura et al.l (|2009l ) for neutrino 
deuteron reactions. However, in this article we first want 
to investigate their possible appearance and their contri- 
bution to the EOS. 

2.3. Model of the HS EOS 

The HS model consists of an ensemble of nucleons 
and nuclei in NSE, whereas interactions of the nucle- 
ons and excluded volume corrections for the nuclei are 
implemented. Note that in the present version of the 
HS EOS tables, the formation of nuclei is neglected for 
T > 20 MeV for simplicity, i.e., matter is assumed to be 
uniform. At such high temperatures there is only a lim- 
ited density region slightly below nuclear matter satura- 
tion density rfi B , where nuclei appear at all. Still it would 
be better to avoid such a hard switch in the EOS, which 
we plan to do in future releases of the EOS tables. In 
the following, we give a brie f summary of the HS model, 
all de tails can be found in iHempel ~fc Schaff ner-Bielichl 

f|2oToT) . 

For the unbound interacting nucleons (neutrons and 
protons) an RMF model is applied. Its Lagrangian is 
based on the exchange of the isoscalar-scalar er, the 
isoscalar-vector w and the isovector-vector p mesons be- 
tween nucleons. In comparison with more sophisticated 
theoretical approaches, based on nucleon-nucleon scat- 
tering data, it has been recognized that nonlinear a and w 
self-coupling terms are necessary to achieve a reasonable 
description of the EOS at high densiti es and the proper- 
ties of nuclei at the same time, see, e.g. JSugahara fc Tokil 
(1994). An alternative approach, which is not used in 
the present investigation, is R MF mo d els w ith density- 
dependent couplings, see, e.g., iTypell ([20051 ). The free 
parameters of the Lagrangians, the meson masses and 
their coupling strengths, have to be determined by fits 
to experimental data. 

We a pply three different RM F parame terizations : 
TM1 by ISugahara fe Tokl (fT99l. TMA by iToki et al.l 
(fl995[) . and FSUgold by iTodd-Rutel fe Piekarewiczl 
(2005). The Lagrangians of TM1 and TMA have the 
same form and include nonlinear terms of the a and w 
mesons. TM1 was developed together with TM2, which 
were fitted to binding energies and charge radii of light 
(TM2) and heavy nuclei (TM1). TMA is based on an in- 
terpolation of these two parameter sets. The coupling pa- 
rameters gi of the set TMA are chosen to be mass-number 
dependent of the form gi — en + bi/A 0A , with a,i and hi 
being constants, to have a good description of nuclei over 
the entire range of mass numbers. For uniform nuclear 
matter the couplings becom e constants and are given by 
dj. T MA was also used in He mpel fc Schaffner-Bielichl 
(2010), where the HS model was introduced. FSUgold 
includes the coupling between the lu and the p meson in 
addition. This leads to a better description of nuclear 
collective modes, the EOS of asymmetric nuclear mat- 
ter and a different densit y dependence of the symmetry 
energy (Piekarewicz 2007) which is very low at large den- 
sities. The coupling constants of FSUgold are fitted to 
binding energies and charge radii of a selection of magic 
nuclei. 

The HS EOS tables take into account the exper imen- 
tal data on nuclear masses from I Audi et al.l (f2003). For 



the masses of the experimentally unknown nuclei differ- 
ent theoretical nuclear structure calculations in form of 
nuclear mass tables are used. The TMA interaction s 
are combined with the mass table of IGeng et ail (pOOl . 
which is also calculated with the RMF model TMA. Thus 
all nuclear interactions are consistent. This mass table 
lists 6969 even-even, even-odd, and odd-odd nuclei, ex- 
tending from ^O to iocjFm from slightly above the pro- 
ton to slightly below the neutron drip line. The nuclear 
binding energies are calculated under consideration of 
axial deformations and the pairing is included with a 
BCS-typc (5-force. For the parameterization TM1 we do 
not have a suitable mass table at hand, thus we cannot 
avoid the minor inconsistency to also use the table of 
IGeng et al.l (|2005| ). which is based on the TMA parame- 
terization. For FSUgold we take a mass table calculated 
by X. Roca-Maza, which was previously applied for the 
outer crust of neutron stars (|Roca-Maza fc Piekarewiczl 
2008). This table contains 1512 even-even nuclei, from 
the proton to the neutron drip, with 14 < A < 348 and 
8 < Z < 100. Odd nuclei are not included in this table. 
The nuclei were calculated only with spherical symmetry 
and the pairing is introduced through a BCS approach 
with constant matrix elements. The constant matrix ele- 
ment for neutrons has been fitted to reproduce the exper- 
imental binding in the tin isotopic chain and the constant 
matrix element for protons to the experimental binding 
in the N — 82 isotonic chain. 

To describe nuclei in the supernova environment, we 
not only need binding energies, but have to account for 
medium and temperature effects. For the screening of 
the Coulomb field of nuclei in the uniform background 
of electrons we use the most basic expression: for each 
nucleus we assume a spherical Wigner-Seitz (WS) cell at 
zero temperature. More elaborated approaches for the 
Coulomb energy of a multi-compo nent plasma at finite 
tempe r ature can, e.g., be fou nd i n iNadvozhin fc Yudin 
(2005). |Potekhin et al.1 (|2009l) . and lPotekhin fc Chabrier 
(2010). However, we leave this for future studies as the 
Coulomb energy becomes only important at low temper- 
atures so that the simplest expression is sufficient for our 
purposes. 

Finite temperature leads to the population of excited 
states of nuclei. Here, we use the temperature d epen- 
dent degeneracy function of iFai fc Randrupl (jl982). It is 
the same analyti c expression as in the original refer ence 
of the HS model (IHempel fc Schaffner-Bieu^ l2010l) and 
used previously by llshizuka et al.l (|2003D . "but now we 
consider only excitation energies below the binding en- 
ergy of the corresponding nucleus, in order to represent 
that the exci ted states still have to be bound (see, e.g., 
IR5pkd (fl98 4) ) . We note that the inclusion of excited 
states up to infinite energies had only a minor influence 
on the composition but would lead to an unphysically 
large contribution of the excited states to the energy den- 
sity and entropy at very high temperatures. 

We describe nuclear matter as a chemical mixture of 
the different nuclear species and nucleons. As we distin- 
guish between nuclei and the surrounding interacting nu- 
cleons we still have to specify how the system is composed 
of the different particles. Our thermodynamic model is 
built on two main assumptions: first, we assume for un- 
bound nucleons that they are not allowed to be situated 
inside of nuclei, whereas nuclei are described as uniform 
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hard spheres at saturation density n B . Second, for nu- 
clei (with mass number A > 2) we assume that they must 
not overlap with any other baryon in the system (nuclei 
or unbound nucleons). Thus, we take the volume which 
is available for the nucleons to be the part of the total 
volume of the system which is not excluded by nuclei. 
This is described by the filling factor of the nucleons 

Z = l-Y,An AiZ /n° B (7) 

A,Z 

(here and in the following, we refer to A > 2). The free 
volume in which a nucleus can move is the total volume 
minus the volume filled by nuclei and nucleons. This is 
incorporated via the free volume fraction 

K = l - n B /n° B , (8) 

with the total baryon number density n B , which includes 
contributions from unbound neutrons and protons: 

n,B =n n + n p + ^2 A tia,z ■ (9) 
A,z 

Based on these two main assumptions, the EOS is 
derived in a consistent way, using the non-relativistic 
Maxwell-Boltzmann description for nuclei and the full 
Fermi-Dir ac integra l s for nucle ons (solved with th e rou- 
tines from lApariciol (|1998h and iGong et~aTl (|2001h ). We 
obtain modifications of all thermodynamic quantities due 
to the exclud ed volume, in agreement with the results of 
lYudinl ()2011[ ) who studied excluded- volume schemes in 
general. Here we give the thermodynamic potential, the 
free energy density /, as an example: 

f = J2fA.z(T,n A ,z)+J2fA:z l 

A,Z A,Z 

-Tj2n A ,zH^ + ^KMF(T,n n /^n p /^ .(10) 

A,Z 

The first term in Eq. (110[) is the summed ideal gas ex- 
pression of nuclei. The Coulomb free energy of nuclei 
appears in addition. The third term is the direct contri- 
bution from the excluded volume. Because of this term, 
as long as nuclei are present, the free energy density goes 
to infinity when approaching saturation density, because 
the free volume of nuclei goes to zero, k — > 0. Thus, 
nuclei will always disappear before saturation density is 
reached. The RMF contribution of the nucleons, /{j MF , 
is weighted with their filling factor £, as the free energy 
is an extensive quantity. If nuclei are absent, £ = 1, and 
we get the unmodified RMF description, as it should be. 
The excluded volume correction for the nuclei represents 
a hard-core repulsion of the nuclei at high densities close 
to saturation density. Instead the modification of the 
free energy of the unbound nucleons is purely geometric 
and just describes that the nucleons fill only a fraction of 
the total volume. In this sense, the two aforementioned 
model assumptions for the excluded volume are essential, 
as they lead to the desired limiting behavior of the EOS. 

2.4. EOS Characteristics and Constraints 

Table [5] lists some characteristic saturation properties 
of uniform bulk nuclear matter for the three different 
RMF parameterizations. We also include results for the 



TABLE 2 

Properties at Saturation Density of the Different EOS 
Under Investigation 





a B 


E 


K 


K< 


J 


L 


EOS 


[fm- 3 ] 


[MeV] 


[MeV] 


[MeV] 


[MeV] 


[MeV] 


TM1 


0.145 


16.3 


281 


-285 


36.9 


111 


TMA 


0.147 


16.0 


318 


-572 


30.7 


90 


FSUgold 


0.148 


16.3 


230 


-524 


32.6 


60 


LS (180) 


0.155 


16.0 


180 


-451 


28.6 


74 


LS (220) 


0.155 


16.0 


220 


-411 


28.6 


74 



Note. Definition of the quantities given in the text. 



LS EOS with the incompressibilities of K = 180 and 
220 MeV. The quantities shown in Table [2] correspond to 
the coefficients of the following power-series expansion 
of the binding energy per baryon at T — around the 
saturation point: 

E(x, (3) = -Eq + ^Kx 2 + ^K'x 3 + ■■■ 

+(3 2 (j+Ilx + ---^ + ■•• , (11) 

with x = n B /n B — 1 denoting the relative deviation 
from the saturation density, and the asymmetry pa- 
rameter /3 which is given by the total proton fraction: 
@ = 1 — 2Y p . Next we discuss the listed nuclear matter 
properties briefly and compare with experimental con- 
straints. The saturation point itself, i.e., the saturation 
density n° B and the corresponding binding energy per nu- 
cleon Eq , can be determined from the analysis of nuclear 
binding energies. All EOS models give reasonable val- 
ues for n B and Eq. The nuclear incompressibility K can 
be studied experimentally, e.g., by exciting the isoscalar 
giant monopole resonanc e (ISGMR) , leadin g to values 
of K = 240 ± 10 MeV (|Piekarewicd [20101) . FSUgold 
and LS (220) agree roughly with this experimental con- 
straint, TM1 is still reasonable, but the value of TMA 
is too large and the value of LS (180) too low. How- 
ever, it is perceived in the literature that the extraction 
of K from ISGMR data is not unambiguous. The exper- 
iments probe nuclear matter slightly below saturation 
density and furthermore the deduced results depend on 
the nuclear interactions (e.g., the behavior of the sym- 
metr y energy) whic h are t a ken for the d a ta ana l ysis (see, 
e.g.. IShlomo et al.l (|2006[ ): IPiekarewicz] (|2004D : ISharmal 
(2009)). One also still has difficulties to explain all avail- 
able IS GMR data consis tently with a single theoretical 
model (Piekarcwicz 201(3). The skewness coefficient K' is 
the third-order term for the expansion of the energy per 
nucleon around n B and is not constrained directly by ex- 
periments so far. Still there exist some estimates: based 
on a correla tion study of Skyrme-Hartree-Fock models, 
[Cher] (120121) deduces K' = -355 ± 95 MeV. Only TM1 
and LS (220) lie within this region. 

Typical values for the symmetry energy coefficient J 
range from 28 to 34 MeV, whereas most studies give 
mean values around 32 MeV, see, e.g., iTsang et al.l 
( 2009) . Thus the symmetry energy of TM1 does not com- 
ply with the experimental results, but all other EOS do. 
We remark that the values of J shown for LS, which we 
calculated with the numerical routines which are avail- 
able online, agree only approximately with the value 
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12 14 16 1£ 

Radius R [km] 

Fig. 1. — Mass-radius relation of cold neutron stars for the EOS 
under investigation. The gray and black horizontal lines show the 
measure ment of the mass of th e pulsar PSR J1614-2230 with its 
ler error IIDe morcst ct al. 2 0101 ). 



given in the original work by lLattimer fc Swestvl (|1991[ ) 
of J = 29.3 MeV. The slope of the symmetry energy, 
L, shows an important co rrelation with the neutr on skin 
thickness of heavy nuclei (]Roca-Maza et ai]|201lD which 
is currently being measured by the PREX experiment at 
JLAB. Interestin gly, L is also correlated with t he radii 
of neutron stars (|Horowitz~fc Pickarcwicz 200 1]). There 
are various experimental probes for L, leading to con- 
straints ranging from 40 to 110 MeV with a common 
center a round 75 MeV, see th e co mprehensive compila - 
tions of ICarbone et al i (|2010T > and iTsang et al.l (|2009j) . 
In some nuclear models a large value of J is compen- 
sated by a large value of L in the fitting procedure of the 
free parameters, to achieve a similar symmetry energy 
slightly below saturation density. This is also seen for 
TM1, which gives too large values of J and L. The other 
investigated EOS give reasonable results for L. 

We want to emphasize the following, simple point 
which will reappear several times during the later discus- 
sion of our results: single parameters do not characterize 
the global behavior of the EOS. To give one example, one 
could think that TM1 is softer than TMA because of the 
lower incompressibility. However, the slight change of 
the saturation point and less negative higher order terms 
such as the skewness coefficient in TM1 lead to the un- 
expected result that the pressure of symmetric nuclear 
matter in TM1 is always larger than in TMA, i.e., TM1 
is stiffcr. Furthermore, for supernovae and neutron stars 
the symmetry energy plays an important role. The ex- 
tremely large symmetry energy of TM1 stiffens its neu- 
tron star EOS additionally. Thus for TM1 one obtains 
a maximum mass of cold neutron stars, which is signifi- 
cantly larger than the one of TMA, as we will see next. 

In 2010, the gravitational mass of the millisecond pul- 
sar PSR J1614-2230 was determined to be M G = 1.97 ± 
0.04M (T) by measuring the Shapiro delay (jDemorest et al.l 
2010). It represents the largest robust maximum mass 
measurement at present and thus the most important 
constraint for the mass-radius relations of neutron stars. 
In Fig. [TJ we show this constraint with its la error to- 



gether with the mass-radius relations of the EOS which 
are investigated in this article. The curves have been 
calculated with the tables which are used in the sim- 
ulations presented below by imposing beta-equilibrium 
without neutrinos at practically zero temperature, i.e., 
T = 0.1 MeV, and solving the Tolman-Oppenheimer- 
Volkoff equations. The values of the corresponding max- 
imum masses of neutron stars are 2.219 M , 2.213 M Q , 
2.022 M , 1.739 M , 1.844 M Q , and 2.056 M Q for 
STOS (TM1), HS (TM1), HS (TMA), HS (FSUgold), 
LS (180), and LS (220), respectively. The mass-radius 
curves of STOS (TM1) and HS (TM1) are almost iden- 
tical, because they are built with the same RMF param- 
eterization TM1. This good agreement shows that there 
are only small differences in the pressure-energy density 
relation between the two EOS models, despite the differ- 
ent description of non-uniform nuclear matter. FSUgold 
and LS (180) fail to fulfill the constraint of PSR J1614- 
2230. Despite this, we include FSUgold in the present 
study, to explore the implications of a very soft RMF 
EOS and because of its good description of low-density 
nuclear matter. We also keep LS (180), because it is 
a standard reference for core-collapse supernova simula- 
tions. The stiff TM1 parameterization lies well above the 
maximum mass constraint, whereas TMA and LS (220) 
a re only slightly abov e . 

iSteiner et al.l (|2010D used observational data of masses 
and radii of seven neutron stars with well-determined 
distances in a Bayesian framework to determine nuclear 
matter parameters and the neutron star EOS. The found 
values for the saturation properties of nuclear matter 
agree remarkably well with the aforementioned experi- 
mental constraints. Furthermore, the authors deduce a 
mass-radius relation which is based on the analysis of the 
observational data and find relatively small neutron star 
radii around 11 to 12 km for Mq = 1.4 M Q , indicating a 
soft behavior of the neutron star EOS around saturation 
density. Only LS (180) is compatible with such small 
radii, and FSUgold and LS (220) come close to it. The 
radii of TMA and TM1 seem to be much too large. 

We conclude that TMA has satisfying nuclear matter 
properties, apart from the too large value of the incom- 
pressibility. For TM1 the incompressibility is barely ac- 
ceptable, but the symmetry energy and its slope are too 
large. On the other hand, these two EOS ar e not compat- 
ible w ith small neutron star radii found by ISteiner et al.l 
(2010). FSUgold has very nice nuclear matter properties 
and a lot of success in nuclear structure, but unfortu- 
nately it gives a too low maximum mass. Also LS (180) 
with an incompressibility of 180 MeV is actually ruled 
out by the observation of PSR J1614-2230, and further- 
more its incompressibility and symmetry energy are too 
low. LS (220) gives satisfactory results. It would be good 
if other new EOS tables became available in the future 
which comply better with all the experimental and ob- 
servational results. The emphasis of the present study is 
a different one. We want to explore a broad range of pos- 
sible EOS, including extreme cases, to better understand 
the impact of the EOS in core-collapse supernovae. 

3. RESULTS: GENERAL FEATURES OF THE HS 

EOS 

In the following subsections, we will compare the 
HS EOS with the parameterization TM1 (HS (TM1)), 
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Fig. 2. — Radial profiles of selected hydrodynamic quantities 
with respect to the encl osed baryon m a ss for the 15 Mq pro- 
genitor model from Wooslcy and Weaver (1995) at the progenitor 
stage, comparing the different EOS of HS (TM1) (red solid lines), 
STOS (TM1) (blue dashed lines) and LS (180) (green dash-dotted 
lines). In panel (d), Y e is the electron fraction, and Yj, the electron 
lepton fraction, "light nuclei" in panel (f ) corresponds to the mass 
fraction of alpha particles for STOS and LS, and the mass fraction 
X a of light nuclei with Z < 5 for HS. "heavy nuclei" is the mass 
fraction of the representative heavy nucleus for STOS and LS, and 
the mass fraction Xa of all nuclei with Z > 6 in HS. Panel (h) 
shows the charge and mass number of the representative heavy nu- 
cleus for STOS and LS, respectively, the average charge and mass 
number of the distribution of heavy nuclei for flS. See the text for 
further details. 



introduced in jj 2 , wit h the standard EOS from 
ILattimer fc Swestvl (|1991[ ) with an incompressibility o f 
180 MeV (LS (180)) as well as H. IShen et all (11998a!) 
(STOS (TM1)). Note that STOS is also based on 
the TM1 parameterization. We examine the collapse, 
bounce and early post-bounce pha ses at the example of 
the 15 Mq progenitor model from IWooslev and Weaver! 
(1995). In this section, we do not consider the LS (220) 
EOS or any of the other RMF parameterizations of the 
HS EOS. Thus in the following discussion sometimes we 
omit the specification "(180)" for the "LS (180)" EOS. 
The same holds for the "HS (TM1)" and "STOS (TM1)" 



EOS, where we sometimes use only "HS" or "STOS" for 
simplicity. 

3.1. Progenitor Stage 

The differences between the EOS can already been 
identified at the progenitor stage, illustrated via ra- 
dial profiles in Fig, p for the 15 M Q progenitor of 
IWooslev and Weave r (1995). We note that we take the 
temperature, density, and electron fraction from the pro- 
genitor as the initial configuration for the core-collapse 
simulation, that is why these quantities are the same 
for all three EOS. For the part of the contracting pro- 
genitor star which is shown in Fig. [5J the baryon den- 
sity ranges from 10 10 to 10 5 g/cm 3 and the entropy per 
baryon from 0.5 to 5 k^, corresponding to temperatures 
of 0.7-0.1 MeV. Matter at the center is already slightly 
asymmetric, whereas the electron fraction ranges from 
Y e ~ 0.42 at the very center to 0.5 at the lowest densi- 
ties shown in the figure. Note that the electron lepton 
fraction Yl = Y Ve + Y e , with the net electron- flavor neu- 
trino fraction Y Ue , is on top of the electron fraction, as 
the abundance of electron neutrinos (and all other neu- 
trinos, too) is still negligible. One sees that the com- 
position is dominated by heavy nuclei, only at a region 
around Mb ~ 1 M Q alpha particles contribute with up 
to 10 % to the mass fraction. 

Next we compare the state of the progenitor described 
by the three EOS in detail. We note that the HS EOS 
should give the most accurate description of the progen- 
itor stage, as it includes shell effects and experimentally 
measured binding energies. Medium modifications of the 
properties of nuclei are negligible at such low densities 
and temperatures encountered here. The fraction of light 
nuclei X a with Z < 5 in the HS EOS is almost identi- 
cal to the alpha-particle fraction of STOS, showing that 
the additional light nuclei are negligible at the progen- 
itor stage. LS gives a reduction of the alpha-particle 
fraction of almost one order of magnitude, which can 
be attributed to a well-known err or in the alpha-particl e 
binding energy in the LS routines (jSwestv fc Mvral2 005) . 
In the reference simulations which we show here, we did 
not correct for this error. Apart from this difference, the 
dependence of the alpha-particle abundance on density 
is very similar in the three EOS. 

The fraction of unbound nucleons is much more model 
dependent, due to the interplay between unbound nucle- 
ons and the formation of heavy nuclei. Regarding heavy 
nuclei in the inner layers, there is a clear tendency of 
STOS to give the largest mass numbers, followed by LS 
and HS with the smallest heavy nuclei. As expected from 
stellar nucleosynthesis, iron-group nuclei are found in HS. 
In contrast, STOS gives nuclei around Z ~ 30. Com- 
pared to the composition of the progenitor (not shown 
in Fig. [2j because only limited information was avail- 
able), all EOS show some significant differences, because 
NSE has actually not been reached completely in the pro- 
genitor calculation, even though temperatures are above 
0.5 MeV. For example, there is still up to 20 % of carbon 
in the very center. Also thermodynamic quantities are af- 
fected from the fact that the progenitor is not in full NSE. 
For example, the central entropy per baryon of the pro- 
genitor model (purple circle in Fig. [He)) is slightly larger 
than in HS, because the average mass number of heavy 
nuclei is lower than the value obtained from NSE. LS 
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Fig . 3. — Radial profiles of the entropy per baryon and the nuclear composition for the 15 Mg progenitor model from Wooslcv et al. 
(2002) at the progenitor stage, with the same notation as in Fig. [2] In addition, we also compare with the progenitor data (black dotted 
lines) given by the stellar evolution calculation. Note that only the average mass number, but not the charge number of heavy nuclei, was 
available in the progenitor data and is depicted in the figure. 



gives a similar but a bit larger entropy than HS. In con- 
trast, the entropy of STOS is significantly lower, which 
w e will explain in the di scussion of the 15 M© progenitor 
of IWooslev etaLl (pOOl ) shown in Fig. El 

We discuss this progenitor too, because the avail- 
able data include entropy profiles and more detailed 
information about the nuclear composition. Further- 
more, this progenitor is clo ser to NSE than the one of 
IWooslev and Weaver! (|1995fl shown before. Regarding 
the composition, similar trends of the three EOS as in 
Fig. [2] are observed. In comparison with the original pro- 
genitor data, there is a perfect agreement with HS for the 
mass fractions of the different particles. However, the 
average mass number of heavy nuclei from the progeni- 
tor calculation is still below the NSE value from HS. HS 
and LS reproduce the entropy profile of the progenitor 
very well, but the entropy per baryon of STOS is almost 
0.3 ks lower. For the shown conditions, the entropy is 
mainly given by heavy nuclei. The entropy of heavy nu- 
clei can be split in two contributions: internal excitations 
and kinetic entropy from the thermal movement of nu- 
clei. All three EOS take into account internal excitations 
of nuclei, but describe it in different ways. Furthermore, 
only HS has a distribution of nuclei, whereas STOS and 
LS are based on the SNA. However, the largely reduced 
entropy of STOS is mainly due to the absence of the ki- 
netic entropy of nuclei. The nuclear kinetic entropy per 
baryon is on the order of 0.2 ks in LS and HS, which 
explains the reduction of the entropy in STOS seen in 
Figs. [2] and |3l This means that STOS assumes always 
a lattice of nuclei, even if a crystal would actually form 
only at lower temperatures. Also the energy density, the 
pressure, and the chemical potentials of STOS lack in 
the kinetic contribution of nuclei in the ideal gas regime 
at low densities and moderate temperatures. This is a 
typical problem of SNA models. LS resolves this prob- 
lem by adding an explicit kinetic contribution of heavy 
nuclei to the free energy. 

3.2. Core Collapse and Bounce 

To follow the evolution of the core-collapse phase, 
in Fig. IH we show radial profiles of selected quanti- 



ties for the three EOS under investigation (LS (180), 
STOS (TM1), HS (TM1)), when the central tempera- 
tures reach ~ 1 MeV. The selected states correspond to 
slightly different times in the simulation around 40 ms 
before bounce, but represent similar stages of the evo- 
lution. For the collapse phase until the stage shown in 
Fig. 01 the evolution of the inner core proceeded homol- 
ogously. No shock has formed and any possible change 
of the entropy per baryon is given only by weak reac- 
tions. As long as the densities are low, neutrinos emitted 
can le ave without interactions. As explained in iBruennl 
(1985) and Licbcndorfcr (2005), during the early stage 
of collapse the entropy generation by electron captures is 
roughly balanced by the entropy carried away by neutri- 
nos. Here, we also find that the entropy profiles remain 
almost unchanged. However, a large difference of the 
entropy per baryon developed at ~ 1.4 M between the 
three EOS, which relates to the different transitions from 
non-NSE to NSE. The huge difference between LS and 
the other two EOS is caused by the different non-NSE 
treatment for LS mentioned in § 12.21 It is to some extent 
artificial and can be safely ignored. Due to the ongo- 
ing compression, the matter in the entire collapsing core 
is heated slightly, which in turn increases the fraction 
of alpha particles and nucleons by almost one order of 
magnitude. By comparing the average mass and charge 
numbers of heavy nuclei, one sees that some electron cap- 
tures have taken place, and that the mass numbers have 
been increased. These processes and the evolution occur 
in an overall similar manner for the three different EOS. 
The differences which were observed for the progenitor 
in Fig. [2l i.e., larger nuclei in STOS and LS, different 
nucleon fractions, and a too low entropy in STOS, can 
still be identified. 

However, some quantities which were identical at the 
progenitor stage developed differently. Let us focus on 
the comparison of HS (TM1) and STOS (TM1) which 
are built from the same nuclear interactions. For ex- 
ample, the electron fraction shows an interesting evolu- 
tion. Up to 0.4 Mq, HS leads to a faster deleptonization 
than STOS. We explain this by the u sed pres c riptio n 
of electron captures on heavy nuclei by iBruen n (1985), 
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Fig. 4. — Radial profiles of selected hydrodynamic quantities 
with respect to the enclosed baryon mass du ring the core-collapse 
phase for the 15 Mq progenitor model from Wooslcv and Weaver 
(1995), when the central temperature T centra i ~ 1 MeV. Shown 
are the HS (TM1) EOS at t ph = -42.9 ms (red solid lines), the 
STOS (TM1) EOS at i pb = -42.1 ms (blue dashed lines), and the 
LS (180) EOS at t ph = -39.1 ms (green dash-dotted lines). The 
notation is the same as in Fig. |21 



where the rate decreases with the excess neutrons of a 
nucleus. Even though (Z) / {A) is roughly the same in 
HS and STOS, because nuclei in STOS are significantly 
larger, the number of excess neutrons is larger, leading 
to smaller electron-capture rates. We leave the ques- 
tion whether this effect would persist in the more elabo- 
rated description s of e lectron capture o n hea vy nuclei of 
iHix et al.1 (|2003| ) and iLanganke et~aTl (120031 ) for future 
studies. Between 0.4 M Q < M B < 1 M Q , the delep- 
tonization was slower in HS than in STOS. This is the 
domain which is dominated by electron captures on free 
protons. Because the free-proton fraction is larger in 
STOS than in HS (see Fig. EJe)) electron captures in 
STOS are faster. Because the electron contribution to 
the EOS is dominating in this region, the smaller elec- 
tron fraction of STOS still has a small effect on hydro- 
dynamic variables like the infall velocity or the baryon 
density, which are slightly increased in STOS. 



Figure [5] shows the evolution of the collapse when the 
central temperatures reach roughly 5 MeV corresponding 
to ~ 0.1 ms before bounce. By comparing the fraction of 
light nuclei X a of HS with the alpha-particle fraction of 
STOS in Figs. |4]and[5] (thick lines in graphs (f)) one real- 
izes that other light nuclei than alpha particles have been 
formed with HS in the center during these two stages of 
the collapsing star. The total fraction of light nuclei in 
the center is almost 10% in Fig. [5j whereas the alpha 
particle fraction has decreased below 1%. If weak pro- 
cesses with the additional light nuclei were taken into 
account, in principle this could also modify the neutrino 
transport and deleptonization. However, the fraction of 
light nuclei other than alpha particles starts to rise only 
shortly before the stage shown in Fig. [5] For example, 
the fraction of 3 H gets larger than the alpha-particle frac- 
tion around 2 ms before bounce. At this stage, neutrinos 
are already trapped in the center. Thus we expect that 
the additional light nuclei have only little influence on 
the deleptonization. 

At the stage of Fig. [5] there is a significant number of 
trapped neutrinos. At the center of the star the elec- 
tron and lepton fractions in HS are still lower than in 
STOS. The change of the electron fraction at the center 
during the stages of Fig. 2] and Fig. O i.e., the num- 
ber of electron captures is very similar for the two EOS: 
AYg HS = -0.1261 and AF e s ™ s = -0.1222. Also the 
change of the central lepton fraction was similar for the 
two EOS, they give a similar amount of deleptonization 
in the center. In the outer layers the differences in the 
lepton fraction profiles are less pronounced and it ap- 
pears that STOS and LS are closer to each other than 
STOS and HS. Interestingly, the electron and lepton frac- 
tion profiles show a more step-like change in HS. These 
steps coincide with larger changes of the average mass 
and charge number of heavy nuclei which are caused by 
neutron magic shell effects, as will be further illustrated 
below. 

In the profiles of the entropy per baryon an increase of 
~ 0.5 Vlb is observed in the center. The entropy creation 
occurs between the two stages of Figs. |4] and \E\ before 
neutrinos become trapped and reach weak equilibrium. 
This happens at the stage where the entropy loss by neu- 
trino transport is smaller than the entropy gain from lo- 
cal weak reactions. If we compare the entropy profiles 
of HS and STOS we see that the difference of the en- 
tropy per baryon for Mb > 0.6 M Q is the same as it 
was at the progenitor stage. Contrary, in Fig. [5] at the 
center, HS and STOS have equal entropies per baryon, 
thus there was a larger entropy production in STOS of 
As stos ^o.65k B thaninHSof As HS -0.47k B . Where 
does th is additional e n tropy come from? Based on the 
work of lBruennl (|1985f ). ILiebendorferl (|2005| ) showed that 
the entropy production can be well approximated by 



Tds = -dY e (A(i — E° sc ) , 
A[i = [i e + tip - [i n , 



(12) 
(13) 



where Ef, sc denotes the typical energy of escaping neutri- 
nos which is on the order of 10 MeV, and A/j, > is the 
average chemical energy which is liberated by an elec- 
tron capture. fj, e , fi p , and fi n are the relativistic chem- 
ical potentials of electrons, protons, and neutrons, i.e., 
they include the corresponding rest masses. Note that 



New equations of state in core-collapse supernova simulations 



11 



neutrinos are emitted at a lower energy than they are 
produced, because they thermalize by inelastic scatter- 
ing before escape. As the change of the electron frac- 
tion, the electron neutrino luminosities and the temper- 
ature evolution are very similar in HS and STOS, we 
find that the difference of the change of the entropy 
per baryon, As STOS — As HS ~ 0.2 he, has to originate 
from a difference A/i STOS — A/i HS on the order of sev- 
eral MeV. This means that electron captures in STOS 
give a larger release of entropy than in HS. Indeed, by 
carefully c omparing the chemical potentials of HS and 
STOS (see lHempel fc Schaffner-Bielichl poToh ) one finds 
that the non-relativistic proton chemical potentials (i.e., 
without the proton rest mass) for the conditions encoun- 
tered here are slightly lower in HS than in STOS. This 
is due to the different description of non-uniform matter 
and nuclei. Also light nuclei other than alpha particles, 
which form during the two stages of the collapsing star 
shown in Figs. 2] and [5] (graphs (f)), give a small mod- 
ification of the chemical potentials in HS. On the other 
hand, STOS has in general larger electron fractions in 
the center, so that fi e will be larger, leading to an in- 
crease of A/x STOS compared to A/x HS . This means that 
electron captures are more energetic, because matter in 
STOS is still more symmetric than in HS. 

Furthermore, the nucleon rest masses employed in the 
EOS have an interesting effect on the entropy production 
as we will discuss now. The EOS model of STOS does not 
include the neutron-to-proton rest-mass difference, but 
only a common nucleon mass of 938 MeV. Contrary, HS 
is built on the measured neutron and proton masses. For 
the conditions encountered here, the nucleon rest masses 
contribute directly to the energy and free energy densities 
and the relativistic chemical potentials. Thus, the rest- 
mass difference is included in A/i HS but not in A/x STOS , 
which alone changes A/i already by Am ~ 1.2935 MeV. 
We conclude, that the appearance of light nuclei, the 
faster deleptonization in the early phases of the collapse 
and realistic nucleon masses lead to less entropy gain for 
electron captures in HS. 

We want to emphasi ze that i n the electron/positron 
capture rates based on IBruennl (|1985l) which determine 
dY e , the nucleon rest- mass difference Am is taken into 
account, independent of the EOS. This is done by using 
the chemical potentials without rest mass and explicitly 
adding Am to the Q-value. Thus, in both simulations 
the energy of the produced neutrinos and the rate of 
deleptonization is determined correctly taking into ac- 
count the rest-mass difference of the nucleons. The miss- 
ing rest-mass difference in STOS leads only to a slight 
overestimation of the entropy production, as discussed 
before. This also explains why the change of the electron 
fraction is so similar for the two EOS, despite the differ- 
ent entropy production. In this context we also want to 
remark that the nucleon rest masses have only a negligi- 
ble effect on the equilibrium (NSE) composition, as long 
as the same rest masses are used for all nucleons, i.e., for 
free nucleons and nucleons bound in nuclei. 

Let us turn to the temperature evolution. In Fig.21 the 
central temperature was equal for STOS and HS. Now 
the temperature of STOS has increased to larger values 
than in HS which is in agreement with the larger entropy 
increase in STOS. However, one also has to note that 
the densities reached in HS are slightly smaller than in 
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Fig. 5. — Same as Fig. [2] but at T central ~ 5 MeV. The states 
correspond to t pb = -0.91 ms for HS (TM1), t pb = -0.84 ms for 
STOS (TM1), and t ph = -1.09 ms for LS (180). 

STOS, so there is less compression heating. In this con- 
text we have to remark that minor differences can also be 
due to the different data sampling, because profiles which 
are shown in the same figure can correspond to slightly 
different evolutionary states. However, we checked that 
the differences which we see and which we discuss here 
are systematic and not caused by the different data sam- 
pling. 

Less than 1 ms later, the core bounce takes place, 
which is depicted in Fig. [51 It is defined as the mo- 
ment when the maximum central density is obtained at 
the end of the iron-core collapse. Regarding the Y e and 
entropy profiles we can identify the same differences be- 
tween HS and STOS as described before. From Fig. [5] 
we can see that the shock is located at a slightly lower 
mass in HS. It is well known that the mass of the in- 
ner core is proportional to Y} (se e iGoldreich fc Weberl 
1980; Mar tmez-Pinedo et al.l l2006h . which explains this 
result to some extent. During the collapse evolution the 
temperature of HS was similar to the one of STOS. Now 
at bounce we see that the core temperature of HS has 
become larger than in STOS. It can also be attributed 
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Fig. 6. — Same as Fig. [4] but at bounce, defined as the moment 
when the maximum central density is obtained at the end of the 
iron-core collapse. 

to the lower electron fraction at the center. As elec- 
trons provide a large contribution to the pressure at sat- 
uration density, the pressure in HS is lower, leading to 
more compression and larger temperatures. We empha- 
size again that STOS and HS are identical above sat- 
uration density. For the very soft non-relativistic LS, 
the same differences to STOS appear but are even more 
pronounced: the shock forms at smaller mass, the infall 
velocities are larger, and the central densities are higher. 
The central conditions at bounce, comparing the three 
EOS, are listed in Table whereas s denotes the en- 
tropy per baryon. We conclude that even though HS 
and STOS are based on the same nucleon interactions, 
the different description of non-uniform nuclear matter 
and nuclei leads to a different evolution of the iron-core 
collapse. HS has a lower central electron fraction, which 
results in a more compact configuration at bounce. 

Regarding the composition there are two main differ- 
ences between HS, and STOS and LS: in the mass shells 
where the shock has ran through in HS there is a signif- 
icant fraction of more than 10% of light nuclei present. 
The alpha-particle fraction of LS and STOS is almost 



TABLE 3 

Selected Central Hydrodynamic Quantities at Bounce for 
the 15 M Models 



EOS 


P 




Ye 


M a 




[10 14 g/cm 3 ] 


[k S ] 




[M~ 


LS (180) 


3.948 


1.4463 


0.2860 


0.5799 


STOS (TM1) 


3.061 


1.1512 


0.3120 


0.6586 


HS (TM1) 


3.135 


1.1484 


0.3063 


0.6258 



Note. 

a Baryon mass enclosed inside the shock. 

two orders of magnitude smaller (see Fig. [6jf)). Obvi- 
ously, such a large abundance of light nuclei influences 
the EOS. The second main difference in the composition 
of the different EOS concerns heavy nuclei. At densi- 
ties larger than 10 14 g/cm 3 there is a small region where 
heavy nuclei appear in HS (see Fig. Hth) between 0.2 
and 0.4 M Q ). This can be seen as the beginning of the 
uniform nuclear matter phase and part of the transition 
to uniform nuclear matter. Instead in STOS and LS 
no heavy nuclei form behind the shock. One observes 
a different behavior of the transition to uniform nuclear 
matter in HS in comparison to LS and STOS. This is also 
the density regime, where the formation of the nuclear 
pasta phases is possible. In HS the pasta phases are not 
taken into account, a Maxwell construction is used, and 
furthermore the transition depends crucially on the im- 
plementation of the excluded-volume effects. However, 
at these high densities neutrinos are trapped. Hence, 
the heavy nuclei at high densities should have a negligi- 
ble effect on the neutrino transport especially since the 
heavy nucl ei lie outside th e approximated Gamow-Teller 
window of lBruennl (|1985f K 

3.3. Post-bounce Evolution 

Figure [7] shows the evolution for the three different 
EOS at 200 ms post bounce. The bounce shock stalls for 
all models and turns into the standing accretion shock, 
which stands around 100 km at 200 ms post bounce (see 
Fig. Eta)). It separates the inner high-density and high- 
temperature core from the outer material which is being 
accreted. Note the entropy differences ahead of the shock 
between LS and STOS/HS, which are due to the afore- 
mentioned different treatment of non-NSE. During the 
long-term post-bounce evolution of several hundreds of 
milliseconds, the standing accretion shocks contract in a 
similar fashion for all EOS under investigation and hence 
the profiles become increasingly similar. The differences 
obtained at bounce due to the different composition and 
different description of nuclei remain in the post-bounce 
evolution. Matter is slightly more neutron-rich for HS 
than for STOS, and HS leads to higher central densities 
and temperatures. In comparison to LS, both HS and 
STOS are significantly stiffer. With LS the shock con- 
tracts much faster and the central density and tempera- 
ture are much higher, as well as Y e is lower (see Fig. Eta), 
(b), (c) and (g)). 

Also the composition still shows some similarities to 
the situation at bounce. The infalling matter is com- 
posed of heavy nuclei, which are systematically smaller 
for HS than in STOS and LS. The accreted matter is 
heated up and dissociates partly into alpha particles. 
When matter reaches the accretion front it encounters 
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strong shock heating, so that nuclei get dissolved and al- 
most only free nucleons remain. Further inside the proto- 
neutron star the compression becomes strong enough 
that the fraction of light nuclei increases again. In gen- 
eral, alpha particles and light nuclei are favored in the 
innermost layers not only because of the higher densi- 
ties, but also because of the lower entropies. In HS 
there is a small region slightly below saturation den- 
sity where even heavy nuclei appear, before they dissolve 
into homogeneous nuclear matter at slightly higher den- 
sities. However, the fraction remains small, less than 
20 % (see Fig. [7][f)). Heavy nuclei cannot be found be- 
low the shock with STOS or LS at the shown time. But 
also with STOS heavy nuclei are present inside the core 
with significant mass fractions but only from shortly af- 
ter bounce until roughly 50 ms post bounce, when they 
disappear again. For STOS on the other hand, alpha 
particles extend down to the center of the proto-neutron 
star in Fig. [7f f ) . In this model, due to the different de- 
scription of the excluded volume effects, alpha particles 
can also exist slightly above saturation density. However, 
more sophisticated studies show th at the Mott transi tion 
actually occurs at lower densities (jTvpel et al.ll20l6T ). 

By looking at the electron fraction in more detail, one 
sees that the envelope of the proto-neutron star below 
the standing bounce shock has a larger electron fraction 
with HS than in STOS. This can be attributed to the 
larger symmetry energy of HS in comparison to STOS, 
due to the inclusion of the additional light nuclei besides 
alpha particles in the HS model. Light nuclei exist as ad- 
ditional degrees of freedom for symmetric nuclear matter, 
which lower the chemical potential of symmetric nuclear 
matter and therefore decrease the Q- values of electron 
captures. As charged-current reactions on light nuclei 
are not taken into account in the present investigation, 
this leads to less electron captures than in STOS. We see 
in Fig. [TJf) that the fraction of light nuclei in HS be- 
hind the shock is roughly two orders of magnitude larger 
than in STOS (and LS). At R ~ 20 km the light nuclei 
fraction increases beyond 10 % in the NSE models, but 
jumps to zero between 10 and 20 km, because in the HS 
EOS nuclei are neglected for T > 20 MeV for simplicity. 

3.4. Nuclear Composition from Collapse to Post 
Bounce 

Figure |8j gives a more detailed view of the contribu- 
tion of light nuclei to the EOS. The mass fractions of 
the most important light nuclei with Z < 3 are plot- 
ted for the simulation with HS (TM1). The black lines 
depict the particles which are also included in LS and 
STOS, and the blue lines depict light nuclei, which were 
usually not considered in core-collapse supernova simu- 
lations until now. We want to emphasize that all mass 
fractions which are shown in this article are a direct out- 
put from the EOS applied in the simulations. In Fig. [5J 
one sees that additional light nuclei apart from alpha 
particles are negligible at the progenitor stage. This ex- 
plains also the good agreement for the fraction of light 
nuclei in Figs. [2] and [3] However, at later times during 
the core collapse where temperature and density increase 
sufficiently, other light nuclei in addition to alpha parti- 
cles appear in non-negligible fractions on the order of a 
few percent and they can become even more abundant 
than alpha particles. In general, the favored appearance 
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Fig. 7. — Radial profiles of selected hydrodynamic quantities 
with respect to the radius f or the 15 Mq progenitor model from 
Wooslcy and Weaver (1995) at 200 ms post bounce. The same 
notation as in Fig. [2] is used. 

of weakly bound light nuclei instead of alpha particles is 
driven by the increased entropy of a system composed 
of smaller fragments, where binding energies play only a 
minor role. How does this work? 

Using Maxwell-Boltzmann distribution functions for 
nuclei, schematically the number density of a nucleus of 
mass A is 

n A cc,4 3 / 2 exp(V^±^ , (14) 

where /i is the nucleon chemical potential without rest 
mass and Ba > the binding energy per nucleon of 
nucleus A. We remark that Ba + /i is always negative 
in our calculations. Thus we do not see any indication 
for Bose-Einstein condensation and the use of Maxwell- 
Boltzmann statistics is justified. For any given finite tem- 
perature, in the limiting case of vanishing density, // will 
be so negative that e A ^l T is completely dominating and 
the binding energies are negligible. Only free neutrons 
and protons will be present. This is an entropy effect, 
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Fig. 8. — Radial profiles of the mass fractions of the most im- 
portant light nuclei with Z < 3 in the HS (TM1) EOS for the 
core-c ollapse simulation of the 15 Mq from Wooslcy and Weaver 
(1995), at the progenitor stage (top), during the collapse (middle) 
and at bounce (bottom). "Li" shows the sum of the mass fractions 
of all lithium isotopes. 

as a result of the minimization of the free energy. With 
increasing density the role of temperature decreases and 
at some point it is favorable to form nuclei, because the 
internal energy contribution to the free energy becomes 
dominating over the entropy contribution. Let us imag- 
ine that all nuclei had the same binding energy per nu- 
cleon B. In this scenario, the average mass number of the 
distribution of nuclei would increase continuously with 
density. Thus naturally there would be a certain density, 
where deuterons gave the main contribution to the com- 
position, followed by alpha particles at larger densities. 
For sufficiently large temperatures T > Ba, so that the 
differences in binding energies per nucleon of light nuclei 
can be seen as small corrections, this scenario gives a vi- 
able explanation for the large mass fractions of weakly 
bound light nuclei observed in our simulations. Contrary, 
for low temperatures T < Ba, the strong binding energy 
of the alpha particle is so dominant that other light nuclei 
will never give the main contribution to the composition. 

It is also instructive to realize that for any given fi- 
nite temperature there is always a certain density where 
the mass fractions of deuterons and alphas are equal. 
This happens approximately at a chemical potential of 
[i ~ Bd — 2B a ~ — 14 MeV, where the difference of the 
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Fig. 9.— Mass fraction distribution of the HS (TM1) EOS for 
three different mass shells, and five different times, ranging from 
the progenitor to 50 ms post bounce. The vertical lines show the 
corresponding average mass number of heavy nuclei. 

total binding energies is compensated by the difference in 
total chemical potentials. At higher densities, the alpha 
fraction is larger, at lower densities the deuteron fraction. 
With increasing temperature this point moves to higher 
densities. Thus only for sufficiently large temperatures 
the fractions of deuterons and alphas will be comparable 
to the mass fraction of unbound nucleons. In conclusion, 
the significant appearance of weakly bound light nuclei is 
driven by an interplay of temperature and density effects. 
Let us come back to the discussion of Fig. [H At 1 ms 
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Fig. 10. — Distribution of t he ma ss fractions of nuclei in the HS (TM1) EOS for the mass shell of Mb = 0.6 Mq and four selected times, 
which are also shown in Fig. [9(a)1 Gray filled squares correspond to nuclei which are included in the HS (TM1) EOS, but which have a 
mass fraction below 10 -10 . The crosses mark the average of the distribution of heavy nuclei and the black squares show stable nuclei as a 
reference. 



before bounce, temperatures above 5 MeV are reached at 
the center (see Fig. [5]) which are sufficiently high so that 
light nuclei can form at large densities. From Fig. [S] one 
sees that tritium ( 3 H), alphas ( 4 He), deuterium ( 2 H), and 
lithium (Li) appear in similar concentrations in the inner- 
most layers of the collapsing core. The largest fractions 
are reached for 3 H because matter becomes increasingly 
neutron rich during the core-collapse phase. We want to 
emphasize that the alpha particle is less abundant than 
deuterium, tritium and the summed lithium isotopes in 
the center. Only in the outer layers the composition is 
dominated by alpha particles. 

At bounce (bottom of Fig. [3]) we see that the formation 
of the shock enhances the appearance of unbound nucle- 
ons due to the increased central densities and tempera- 
tures. At this stage, the matter behind the shock, which 
is located at ~ 0.6 M Q , is heated to roughly 14 MeV (see 
Fig. [6]). Almost all heavy nuclei get dissociated, but the 
fractions of light nuclei are increased. Note that the den- 



sities are very large ~ 10 14 g/cm 3 . At low densities there 
would be no light nuclei for such large temperatures any 
more, but here we have a competition between density 
and temperature effects, as described before. At the very 
center of the collapsing stellar core, at densities above 
saturation density, light nuclei disappear completely and 
uniform nucleon matter is reached. 

Regarding the contribution of heavy nuclei, with the 
HS EOS we can study the impact of the commonly used 
SNA. Figure [9] shows the mass distributions of nuclei 
with the HS (TM1) EOS for the (baryonic) mass shells 
of 0.6, 0.8, and 1.0 M at different times. Fig. [10] de- 
picts the distributions in form of nuclear charts only for 
the 0.6 M Q mass shell. We want to concentrate on this 
mass shell hrst. At t he pr ogenitor stage (Fig. [TUt a) and 
red solid lines in Fig. 9(a) I we clearly see that iron-group 
nuclei are most abundant, already mentioned in the dis- 
cussion of Fig. [2j The distribution is still close to the 
valley of stability. During the ongoing contraction in 
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Fig. 11. — Radial profiles of the baryon density (solid lines) and the temperature (dashed lines) in the top panels and the corresponding 
nuclear composition regarding light nuclei (bottom panels) at two selected post-bounce times for the 15 Mq progenitor model from 
Wooslcy and Weaver (1995) using the HS (TM1) EOS. The green vertical lines denote the position of the neutrinospheres for f e (solid 
lines) and for P e (dashed lines). 



the core collapse, electron captures take place and the 
formation of heavier and more neutron-rich nuclei is fa- 
vored. In Fig. fTOT b). we see the importance of closed 
neutron shells, i.e. nuclei with N = 50 and TV = 82 
are most abundant. These nuclei c orresp ond to the two 
peaks of the blue dotted lines in Fig. 9(a) Shortly before 
bounce, the distribution of heavy nuclei becomes very 
broad and goes beyond A = 300, which is nicely illus- 
trated in Fig. riffl e) . Note the different scale in this sub- 
figure. Very exotic nuclei are found and the distribution 
is far away from the valley of stability. The most abun- 
dant nuclei are at the well-known neutron-magic number 
N = 82, but also other neutron-magic numbers (50 and 
126) are strongly populated, visible by the additional 
peaks and bumps i n the distribution shown by the green 
dashed lines in Fig. 9(a) Obviously, such multiple peaks 
and/or bimodal distributions cannot be captured by an 
average heavy nucleus, depicted by vertical lines in Fig. [2] 
and crosses in Fig. [TO] 

Another interesting aspect of heavy nuclei is the largest 
mass numbers which are reached in the simulations. 
iFurusawa et al.l (|2011l ) pointed out that they obtain nu- 
clei with mass number up to 2000 in their EOS and even 
larger nuclei are present in the STOS EOS. Superheavy 
nuclei with A up to 3000 are also fou nd in the Hartree 
calculations by G. iShen et all ()2011bO . Contrary, in HS 
it is impossible to obtain such superheavy nuclei, due 
to the use of nuclear mass tables based on experimen- 
tal data and nuclear structure calculations. At the times 
shown in all the previous plots the maximum average 
mass numbers are around 250, whereas all EOS lie in 



a similar range. However, in Fig. HOf c) one sees that 
the distribution of heavy nuclei in HS reaches the border 
of the mass table. The overall maximum mass numbers 
which we get during the entire simulation are A ~ 900 
for STOS, but only A ~ 300 for HS (TM1) and A ~ 200 
for LS. For all three EOS these overall maximum mass 
numbers are found at similar conditions, at ~ 0.5 ms be- 
fore bounce and Mb ~ 0.1 Mq. From this point until 
bounce the mass numbers decrease, see Fig. [6l In the 
early post-bounce evolution of STOS there is a period 
where heavy nuclei appear again in the core at densities 
around 10 14 g/cm 3 , i.e., in the transition region to uni- 
form nuclear matter. At 50 ms post bounce these heavy 
nuclei have disappeared again, therefore they cannot be 
seen in the previous figures. This appearance of heavy 
nuclei is similar like the peak in the mass fraction of 
heavy nuclei with HS below the shock in Figs. [5] and [71 
but for HS it persists for a much longer time. Further- 
more, in STOS these nuclei are superheavies with mass 
numbers up to 500. In general we find that neutrinos 
are trapped under the conditions where superheavies ap- 
pear, and thus they do not have an effect on the neutrino 
transport. Therefore we do not expect that the limita- 
tion of HS to nuclei with mass numbers up to ~ 300 is 
crucial for the supernova dynamics. 

We turn back to the evolution of the nuclear distri- 
bution of the mass shell shown in Fig. [TUl At bounce, 
the shock forms very close to 0.6 Mq. Almost imme- 
diately after the stage depicted in Fig. [TOT c). the shock 
runs through the selected mass shell and most of the 
heavy nuclei get dissociated (Fig, flffld)). A distribution 
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of light nuclei forms, with mass fractio ns wh ich decrease 
exponentially with increasing A (Fig. 9(a) I. This form 
of the distributions at 1 and 50 ms post bounce indicates 
that the binding energies are not very important any 
more, due to the high temperatures. When the shock 
propagates outward and dissociates the infalling matter, 
the distributions change from the very broad shape with 
peaks due to binding energy and shell effects to a sim- 
ple exponential. The change of the shape of the dis- 
tributions can also not be captured by the average of 
the heavy nuc lei. F or the higher mass shells shown in 
Figs. 9(b) and 9(c) similar features are observed but at 
later times. Furthermore, the mass numbers reached are 
significantly smaller for the outer mass shells, because 
matter is much less compressed. Animated results of the 
nuclear distributions of different mass shells can be found 
at the aforementioned HS EOS Web site. 

Also in the post-bounce evolution light nuclei give an 
important contr ibution to the co mposition, as can be 
seen from Figs. 11(a) and 11(b) Again, the standard 



particles which are usually considered in a core-collapse 
supernova are shown in black, and the additional new 
light nuclei in blue. One can identify three different re- 
gions: the core of the newly born proto-neutron star 
at the center up to R ~ 10 km, the envelope of the 
proto-neutron star up to the standing accretion shock at 
R ~ 200 km, and the infalling matter in the outer lay- 
ers. The core is composed of mainly free nucleons and 
some light nuclei which disappear at densities larger than 
saturation density. The accreted matter is a mixture of 
mainly heavy nuclei, alpha particles and free protons at 
low abundances. However, behind the standing accre- 
tion shock the sudden increase of temperature dissoci- 
ates heavy nuclei into smaller fragments, a mixture of 
nucleons and light nuclei appears. Note that 3 H is the 
most abundant light nucleus in the proto-neutron star, 
whereas 2 H dominates the shock heated low-density mat- 
ter which is more symmetric. The 2 H fraction is usually 
at least one order of magnitude larger than the alpha- 
particle fraction. 

3.5. Impact on Neutrino Heating and Cooling 

The presence of light nuclei can modify neutrino heat- 
ing and cooling. In the present article, we do not in- 
clude any inelastic weak processes on light nuclei, but 
the detailed knowledge of their abundances gives already 
first insight into the possible impact they might have 
on neutrino heating and cooling. The heating region is 
located between the neutrinospheres and the standing 
bounce shock, whereas cooling occurs around the neutri- 
nospheres. To identify th e abun dan ces of light nuclei in 
the two regions, in Figs. |ll(a) and 11(b) we also show 
the neutrinospheres (see the vertical lines, v e : solid lines, 
v e : dashed lines). During the earlier post-bounce phase 
(see bottom panel of Fig. [5J, 2 H and 3 H are as abun- 
dant as protons right behind the still expanding bounce 
shock. There, the inclusion of inelastic weak processes 
with light nuclei, i.e., inelastic scattering and most dom- 
inantly charged-current reactions, may modi fy the neu- 
trino heating. However, already at 50 ms (Fig. 11(a) ) the 
abundances directly below the shock have decreased and 
the neutrino flux is already quite diluted geometrically 
at this distance close to the shock. The most efficient 
heating occurs at smaller radii, where there are much 
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Fig. 12. — Heating-cooling rates, i.e., the net energy deposition 
by neutrinos for the 15 Mq progenitor models at 200 ms post 
bounce, comparing the different EOS under investigation. The 
vertical lines show the neutrinospheres. 



less light nuclei than free protons. 

The abundance of light nuclei in the heating region de- 
creases further on the long term. The slight co mpres sion 
of the surface of the PNS which is seen in Figs. 11(a) and 
11(b) decreases the density close to the shock, leading to 
a dissolution of light nuclei into nucleons. At 200 ms post 
bounce, above the neutrinospheres only few light nuclei 
are found in comparison to the free nucleons. Hence a 
very strong impact on neutrino heating is not expected 
during the later post-bounce phase s. Nevertheless it is 
interesting to note that in Fig. 11(b) the largest difference 



between the alpha-particle fraction and the fractions of 
additional light nuclei is observed in the heating region 
between the neutrinospheres and the standing bounce 
shock. There, the fraction of alpha-particles is up to 
three orders of magnitude smaller than the fraction of 
deuterium. The largest contributions of light nuclei in 
the post-bounce evolution are found below the neutri- 
nospheres. There, e.g., deuterium is still as abundant as 
protons. It may modify cooling and could have an effect 
on the explosion dynamics. For example, less cooling has 
a similar potential than more heating. Additionally, an 
impact on cooling is interesting because it modulates the 
neutrino signal. Furthermore, it may be of particular in- 
terest for studies that investigate the long-term evolution 
of proto-neutron stars, during which the neutrinospheres 
move continuously to higher densities where light nuclei 
are more abundant. The conclusions about the possi- 
ble impact of light nuclei on the supernova dynamics are 
only qualitatively and should be investigated in more de- 
tail (including inelastic weak processes with light nuclei) 
in future studies. 

Apart from the composition also the EOS affects the 
heating and cooling via the compactness of the proto- 
neutron star. This is shown in Fig. [12] where we plot 
the heating-cooling rates for v e (left panel) and v e (right 
panel) with respect to the baryon density, at 200 ms post 
bounce. It becomes clear that the largest heating be- 
low the standing bounce shock (which is located around 
p ~ 10 8 g/cm 3 ) but also the largest cooling rates at the 
neutrinospheres are obtained for the soft LS (180) EOS. 
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Fig. 13. — Evolution of the neutrino luminosities and root-mean-square energies for the 15 Mq progenitor model from Wooslcy and Weaver 
ll995h , comparing the HS (TM1) EOS (red solid lines), the STOS (TM1) EOS (blue dashed lines), and the LS EOS (dash-dotted green 
lines). The observables are sampled in a co- moving reference frame at a distance of 500 km. Note the two different scales in the subfigure 
for the electron neutrino luminosities. 



Although STOS (TM1) and HS (TM1) are both based 
on the same RMF parameterization TM1, the slightly 
more compact proto-neutron star of HS leads to larger 
heating and cooling rates for HS. 

3.6. Neutrino Signal 



Figures 11(a) and ll(b)| show that neutrons and pro- 
tons are the most abundant particles at the neutri- 
nospheres, but also 2 H and 3 H appear at non-negligible 
abundances of lO -1 to 10~ 2 . Consequently, the inclu- 
sion of the weak processes with light nuclei which are 
not incorporated in the present study may modify the 
neutrino signal. Here, we will identify the impact of the 
three EOS on the neutrino signal. Several aspects re- 
lated to the comparison of the neutrin o signal of STOS 
and L S ha ve already been disc ussed in iSumiyoshi et all 
(2008) and lFischer etall (pOOfll ) for various massive iron- 
core progenitor models from different groups. Here we 
will add the differences and similarities between HS and 
STOS as well as LS to the discussion. 

Figure [T3l depicts the neutrino signal, i.e., the evolu- 
tion of luminosities and mean energies, for the 15 M Q 
progenitor comparing the three reference EOS under in- 
vestigation. Note the two different scales in the sub- 
figu re for t he electron neutrino luminosity (top panel of 
Fig. 13(a) ). Let us first look at the collapse phase before 
bounce during which only v e are produced via electron 
captures at heavy nuclei and free protons. The earliest 



time which is shown in Fig. 1131 corresponds to a stage be- 
tween Fig. S] (t pb ~ —40 ms) and Fig. [5] (t ph ~ — 1 ms). 



If we look at the top panel of Fig. 13(b) we see that the 
mean energies of electron neutrinos are about 1.5 MeV 
smaller comparing HS with LS and about 2 MeV smaller 
comparing HS with STOS. This is related to the con- 
ditions at decoupling at the neutrinospheres R v during 
the core-collapse phase, where we find the following or- 
dering: pg s > p| TOS > p^S > t hs > t ls > T s^o S; 

and iZHS < flsjos „ flLS 

Since low-energy neutrinos 
decouple generally at high densities and temperatures 
with small neutrinosphere radii, the following hierarchy 
holds for the mean energies (E)^ < (E)^ < (2?)f£ • 
The electron neutrino luminosities which are also shown 
in Fig. 
the col 



13(a) are determined by mass accretion. During 



apse phase the differences in the mass accretion 
rates are small for the different EOS and the ordering 
of the neutrinospheres leads to the following hierarchy 
for the luminosities L^ s < L^ TOS < L]f. During the 
first 10 ms after bounce, when the shock runs through 
the neutrinospheres, the ^ e -deleptonization burst is re- 
leased. We find that the luminosities and mean energies 
of LS are lowest with maximum L Ve ~ 4 x 10 53 erg/s, 
whereas HS and STOS are very similar with maximum 
L„ o ~ 4.5 x 10 53 erg/s. 

During the later post-bounce evolution, v e and v e 
are produced mainly via electron captures at protons 
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and positron captures at neutrons. Their luminosi- 
ties are given by mass accretion at the corresponding 
neutrinospheres. The following hierarchy holds for the 
mass accretion rates at the neutrinospheres to^ s _ > 

m^ s > ?rjfj™ s . Note also the ordering of the neutri- 

nosphere radii, i?^ s < i?^ s < Rf TOS . The different mass 
accretion rates are related to the different post-bounce 
conditions obtained for the different EOS. LS leads to the 
most compact configuration where the standing bounce 
shock contracts fastest which in turn leads to the highest 
mass accretion rate. The opposite holds for STOS, and 
HS lies between LS and STOS. Note that the small dif- 
ferences between HS and STOS are due to the different 
description of light and heavy nuclei at low and inter- 
mediate densities, which leads to slightly lower Y e and 
a slightly more compact configuration for HS. Hence, 
^vjue > ^V«%e — throughout the entire post- 



TABLE 4 

Selected Central Hydrodynamic Quantities at Bounce for 
the 40 M Models 



bounce phase shown in Fig. 13(a) up to 650 ms. On the 
other hand, Z/W T are produced by pair processes and they 
interact only via neutral-current reactions. They de- 
couple at generally higher densities. Consequently their 
neutrinospheres are deeper inside the proto-neutron star 
than the electron flavor neutrinospheres. The v^/ T lumi- 
nosity can be expressed in terms of diffusion which scales 
as T 4 \]i v ^ . Here, higher temperatures at the neu- 

trinospheres are compensated by smaller neutrinosphere 
radii and hence the luminosities have similar values for 
all EOS under invest igation during the post-bounce evo- 
lution shown in Fig. 13(a) up to 650 ms. 

The ordering of the mean energies of the different 
EOS is the same as for the luminosities. It reflects 
the ordering of the temperatures at the neutrinospheres. 
However, mean energies of (/i/r)-neutrinos for LS are 
significantly higher than for HS and STOS. Note that 
(pi/r)-neutrinos are sensitive to differences at high den- 
sities, where the EOS under investigation lead to differ- 
ent proto-neutron star contraction behaviors. The sig- 
nificantly faster proto-neutron star contraction for LS, 
in comparison to HS and STOS, leads to the h igher 
(pi/r)-neutrinos mean energies shown in Fig. 13(b)| The 
contraction behavior for HS and STOS is more similar 
at high densities, leading to almost equal mean (p/r)- 
neutrino energies. 

4. RESULTS: DIFFERENT NUCLEAR 
INTERACTIONS 

The EOS imprint on the neutrino signal for the 15 
Mq progenitor is rather weak. However, if we go to 
the regime of failed supernovae of very massive progeni- 
tors, the influence of the EOS becomes more pronounced, 
as was also found i n prev io us EOS studi e s, like , e.g., 
bv iSumivoshi et all (120071) . iFischer et all (|2009l ). and 
IQ'Connor fc Otti(|2011h~ In order to illustrate and under- 
stand the effect of the high-density EOS on the emitted 
neutrino signal from core-collapse supernovae, we simu- 
late the collapse of the ma ssive 40 M Q progenitor from 
iWooslev and Weaverl ()1995l ) which eventually ends up in 
a black hole. For such a massive progenitor, very high 
densities and temperatures are reached, where different 
EOS show significant differences. We want to empha- 
size that this is the first time that the non-relativistic LS 
EOS are compared with several different RMF EOS, giv- 



EOS 


P 




Y e 




[10 14 g/cm 3 ] 


[k B ] 




LS (180) 


3.855 


1.487 


0.2844 


LS (220) 


3.674 


1.475 


0.2826 


HS (FSUgold) 


3.151 


1.529 


0.2815 


HS (TMA) 


2.878 


1.545 


0.2811 


STOS (TM1) 


2.690 


1.505 


0.2915 



ing a more comprehensive view of EOS effects. On the 
other hand, we do not explore the important sensitivity 
on the progenitor star w hich was thoroughly studied by 
IQ'Connor fc"Ottl (IMlTl) . 

From the relativistic EOS models we apply 
STOS (TM1) and the two new EOS tables HS (TMA) 
and HS (FSUgold). Because HS (TM1) and 
STOS (TM1) are identical at high densities, we do 
not include HS (TM1) in the comparison any more. 
The two new EOS tables HS (TMA) and HS (FSUgold) 
have been introduced in §2.31 and §2.4| and selected 
EOS properties are listed in Table El We consider two 
different versions of the non-relativistic LS EOS in the 
comparison, LS (180) and LS (220), characterized by 
different values of the incompressibility (see Table [2]). 
However, the simulation with LS (220) is performed 
with a lower resolution where we use only four neutrino 
propagation angles, because for this model we restrict 
the discussion on hydrodynamical aspects. 

The most important hydrodynamical quantities at the 
core bounce of the five different EOS are given in Table|4j 
If one compares these values with the EOS properties of 
Table [21 one sees that LS (180) with the lowest incom- 
pressibility has the highest central density at bounce, 
followed by LS (220) and then FSUgold which has an 
incompressibility of 230 MeV. For TM1 and TMA the 
correlation between the density at bounce and the in- 
compressibility is inverted. TM1 has the lowest density, 
even though the incompressibility of TMA is the largest. 
This has two reasons: here, we have asymmetric nuclear 
matter and the large symmetry energy of TM1 stiffens 
the EOS. Furthermore, higher order terms for the ex- 
pansion around saturation density for symmetric matter 
(e.g., the skewness) are larger in TM1, so that TM1 be- 
haves stiffer, as also noted before in § 12.41 This example 
illustrates that it is generally not possible to explain the 
global properties of an EOS with single saturation pa- 
rameters, especially when going to high densities and/or 
asymmetry. 

The most interesting differences occur in the late post- 
bounce phase, shortly before collapse to a black hole. 
Here, one could expect that the maximum mass of a cold 
neutron star is correlated with the time until black hole 
formation and the correspondin g neutrino signal, as wa s 
also found in p reviou s stud ies of ISumivoshi et all (120071) . 
IFischer et all IpOOflh . and IQ'Connor fc Otti (|2011D . If 
this expectation was true, we should find the follow- 
ing sequence for the time until black hole formation: 
HS (FSUgold), LS (180), HS (TMA), LS (220), and 
STOS (TM1), in increasing order. Next we investigate 
the neutrino signal and try to identify such correlations. 

The neutrino signals for the 40 M Q model using differ- 
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Fig. 14. — Evolution of the neutrino luminosities and mean energies for the 40 Mq progenitor model from Wooslcy and Weaver (1995), 
comparing the HS (TMA) EOS (red solid lines), HS (FSUgold) EOS (black dashed lines), the STOS (TM1) EOS (blue dashed lines), and 
the LS (180) EOS (dash-dotted green lines). The observables are sampled in a co-moving reference frame at a distance of 500 km. Note 
the two different scales in the subfigure for the electron neutrino luminosities. 



ent EOS are shown in Fig. [HI similar to Fig.rjJ LS (220) 
is not included in this figure, due to the reduced neutrino 
angular resolution in its si mulation . We remark that the 
mean energies shown in Fig. 14(b) are always continuous. 
However, in the first milliseconds after bounce, they can 
change so quickly that this cannot be resolved in the 
scaling of the figure. The differences obtained during the 
core-collapse phase before bounce for the 40 M models 
are equivalent to what we have already discussed in the 
previous section at the example of the 15 M models. For 
the two HS models the neutrino light curves and evolu- 
tion of the mean energies lie almost on top of each other, 
independent from the parameterization, until a certain 
level of compactness is reached. This also holds for the 
HS (TM1) EOS which is not shown here. It is very inter- 
esting that all HS EOS give a very similar neutrino signal 
even during the early post-bounce phase. The neutrino 
signal is generated at densities below saturation density, 
and even though the nascent proto-neutron star may be 
different for the different EOS, the imprint of the high- 
density phase is rather weak. Obviously, at low densities 
the interactions of unbound nucleons are not so impor- 
tant, whereas the description of nuclei and the degrees 
of freedom in the nuclear composition are crucial, which 
explains why the neutrino signal evolves so similar for 
the different HS models. 



However, the moment of black hole formation is set 
by the nuclear interactions in the high-density region of 
the EOS. Figure [T5l depicts the evolution of the central 
densities for the five different EOS. After bounce, the 
curves first show a linear behavior, until the newly born 
proto-neutron star becomes gravitationally unstable and 
collapses to a black hole. LS (180) shows a very fast com- 
pression in the center, followed by LS (220), then FSUg- 
old, then TMA and finally TM1 with the slowest com- 
pression behavior. By approaching the moment of black 
hole formation a significant rise of the v^/ T luminosities 
is observed. The different EOS have different timescales 
until black hole formation, and thus pronounced differ- 
ences occur for the v^/t luminosities and mean energies, 
making /Lt/r-neutrinos the best messengers of the EOS. 
Contrary, there is only little separation for the evolution 
of v e and v e mean-energies for the four EOS, and also the 
luminosities remain more similar. Further details about 
the rise of the v^j T luminosities and its connection to 
th e contraction of th e proto -neu tron star can be fou nd 
in lLiebendorfer et all (|2004[ ) and lFischer et all (|2009t ). 

Table [5] lists selected hydrodynamic quantities for the 
last stable configuration before the beginning of black 
hole formation. Here this is defined as the moment when 
infall velocities on the order of 1000 km/s are obtained 
in the PNS interior. We study the possible correlation 
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Fig. 15. — Evolution of the centr al density for the col l apse o f 
the 40 progenitor model from Wooslcv and Weaver (1995), 

comparing the different EOS under investigation. 

between the maximum mass of a cold neutron star and 
short/long accretion time until black hole formation by 
comparing i p b and Mq from Table [5] with the maxi- 
mum masses of cold neutron stars given in § 12.41 The 
results are shown in Fig. I17[ where the black crosses 
depict the time until beginning of black hole formation 
as a function of the maximum gravitational masses of 
cold neutron stars. The maximum gravitational masses 
from Table [5] which were found in the simulations are 
shown by the blue squares. Compared to the cold con- 
figurations they are significantly increased for all EOS, 
but not equally strong. Surprisingly, although LS (180) 
has a larger maximum mass for a cold neutron star than 
FSUgold, and the maximum mass of LS (220) is larger 
than for FSUgold and TMA, the two LS EOS lead to the 
earliest beginning of black hole formation. In the super- 
nova environment they behave less stable than the RMF 
EOS and even FSUgold gives a later black hole forma- 
tion then both of them. FSUgold is followed by TMA 
and then comes TM1 with the latest collapse. Conse- 
quently, the expected ordering with the maximum mass 
of cold neutron stars established in the literature so far is 
not valid any more for our extended set of EOS including 
the new TMA and FSUgold EOS. The hierarchy between 
the maximum mass of cold neutron stars and short/long 
accretion times before black hole formation holds only 
separately within the class of RMF models or within the 
class of the non-relativistic LS EOS. 

We want to examine this result further. In Tabled the 
time until black hole formation and the enclosed baryon 
mass have the same ordering, which shows that the ac- 
cretion rate is affected only little by the EOS. On the 
other hand the central densities at the onset of collapse 
are in most cases lower for the EOS which have a large 
mass at a later collapse. It illustrates the stiffness of the 
EOS for the conditions encountered here. Interestingly, 
for the central temperatures at the onset of collapse to a 
black hole, which are listed in Table El we find roughly 
50 MeV for all three RMF EOS, but only 30 MeV for the 
two non-relativistic LS EOS. 

The hydrodynamic state corresponding to the last sta- 
ble configuration is further illustrated in Fig. 1161 The 



TABLE 5 

Selected Quantities at the Onset of Collapse to a Black 

Hole 



EOS 


-t a 








M G Z 




M 


[10 15 g/cm 3 ] 


[MeV] 


[M ] 


[Mq] 


LS (180) 


0.415 


1.292 


29.978 


2.227 


2.133 


LS (220) 


0.521 


1.324 


31.446 


2.350 


2.233 


HS (FSUgold) 


0.571 


1.058 


48.104 


2.465 


2.341 


HS (TMA) 


0.737 


0.943 


46.708 


2.626 


2.466 


STOS (TM1) 


1.028 


0.769 


49.705 


2.864 


2.652 



Notes: 

a Time post bounce. 

b Baryon density in the center. 

c Temperature in the center. 

d Baryon mass enclosed inside the shock. 

e Gravitational mass enclosed inside the shock. 



similar behavior of the central temperatures of the RMF 
EOS and the different behavior of LS can be seen clearly. 
The peak temperatures around M — 0.8 M Q are dif- 
ferent for all five models. Interestingly, they show the 
same ordering like the central densities at core bounce 
in Table H namely, LS (180), LS (220), HS (FSUgold), 
HS (TMA), STOS (TM1) in decreasing order. This is 
also the sequence for the time until black hole forma- 
tion. Note that the same ordering is also seen in the 
luminosities and mean-energies of /z/T-neutrinos, which 
are most sensitive to temperature effects. The entropy 
profiles look more similar than the temperature profiles, 
whereas the two LS EOS lead to slightly higher entropies 
for Mb > 0.5 Mq. In Fig.fTolone can also identify the dif- 
ferent central densities and different masses of the proto- 
neutron stars mentioned before. The electron fraction 
profiles show only small differences for the investigated 
EOS. The neutrino abundances are rather small and are 
also very similar for the different EOS. 

We come to the unsatisfying conclusion that nei- 
ther the incompressibility (comparing HS (TMA) and 
STOS (TM1)) nor the maximum masses of cold neutron 
stars can directly be related to the time until black hole 
formation. Single saturation properties are not sufficient 
to describe the behavior of the EOS in the simulation, 
where the different EOS evolve to different thermody- 
namic states. The maximum mass of cold neutron stars 
is not very meaningful, because the found states are very 
different from the cold configurations and because the 
strength of the temperature effects is model dependent 
as we will discuss in the next paragraph. 

The profiles at the beginning of black hole formation 
from the simulations shown in Fig.[16]can roughly be ap- 
proximated by a constant entropy per baryon of s = 4 k^ 
and electrons in beta-equilibrium, where the contribu- 
tion of neutrinos is negligible. We can compare the re- 
sults from the simulations with hydrostatic configura- 
tions of the Tolman-Oppenheimer-Volkoff equations by 
using these approximations. The corresponding maxi- 
mum masses are plotted in Fig. [T7] versus the time until 
black hole formation, tmi, from Table 03 by red circles. 
If one compares with the gravitational masses from the 
simulations which are depicted by blue squares, one sees 
that there is a rough agreement between "s = 4" and 
"sim" . Furthermore, the time until black hole formation 
is monotonically increasing with M max (s = 4), similar 
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hole o f the 40 Mq progenitor model from Wooslcy and Weaver 
(1995), comparing the different EOS under investigation 



as for M max (sim). This shows that the above approxi- 
mations (s — 4 ks, beta-equilibrium without neutrinos) 
work sufficiently well to describe the states at the be- 
ginning of black hole formation from the simulations. 
In conclusion, the significant increase of the maximum 
mass in the simulation and the weak correlation between 
-^max(cold) and £bh is explained mainly as a tempera- 
ture effect. 

Very often for hydrostatic configurations of proto- 
neutron stars a constant entropy per baryon of s = 2 k^ 
and electron lepton fraction Yl — 0.4 is assumed. These 
conditions are clearly different to the profiles shown in 
Fig. [16] and are more appropriate for proto-neutron stars 
of less massive progenitors. We remark that such hydro- 
static configurations with s = 2 ks and Yl = 0.4 are in- 
deed not suitable to explain the observed behavior. Their 
maximum masses are 0.3-0.5 Mq lower than observed in 
the simulations, and do not show a monotonic correlation 
with £bh- A large entropy of s = 4 k# is the important 
point for the explanation, whereas the neutrino fraction 



is less significant. 

It is interesting that the non-relativistic LS EOS show 
less stiffening at finite entropy than the relativistic EOS. 
The different temperature profiles discussed before also 
indicate a different temperature dependency of the LS 
EOS. It is difficult to pin down whether this is an ar- 
tifact of the non-relativistic dispersion relation or just a 
result of the chosen parameterization of the nuclear inter- 
actions. For example, the temperature-entropy relation 
depends on the effective nucleon mass m* which appears 
in the dispersion relation. In LS the effective mass is set 
equal to the vacuum mass, m*/m — 1, whereas m*/m 
goes down to 0.2 in the core of the proto-neutron stars 
for the RMF models. The neglect of the effective mass 
reduction at finite density in LS leads to a suppression 
of temperature effects. In addition to the effective mass, 
also the electron fraction affects the temperature-entropy 
relation, because the electron contribution depends dif- 
ferently on temperature than the nucleon kinetic contri- 
bution. The electron fraction on the other hand is set 
by the symmetry energy. These effects, which are inde- 
pendent from the used energy-momentum relation, can 
be identified by comparing different RMF models or the 
two LS EOS. For example, LS (220) shows less depen- 
dence on entropy than LS ( 180). The same result was 
found by lO'Connor fc~Ottl 0W i. It could be due to 
the dominance of the LS (220) EOS by interactions com- 
pared to the kinetic contribution. For the RMF models 
we see that FSUgold has the largest increase of the max- 
imum mass at constant entropy. This can be related to 
its softer density dependence of the symmetry energy at 
high densities. 

Independent of the detailed origin, with these results 
we have found a quantity which is directly correlated 
with the time until black hole formation: the maxi- 
mum mass of proto-neutron stars, approximated by the 
"s = 4" configuration. The time until black hole for- 
mation is monotonically increasing with the maximum 
masses of "s = 4" proto-neutron stars with an approxi- 
mately linear dependency. If for a known progenitor one 
could detect and follow the neutrino signal of the core- 
collapse supernova until black hole formation where the 
neutrino signal ceases, the correlation of Fig. [T7] would 
give information about the maximum masses of such 
"s = 4" configurations. This would represent a signif- 
icant constraint for the nuclear EOS at finite entropy, 
complementary to observations of cold neutron stars. 

The mass measure ment of PSR J1614-2230 by 
iDemorest et~aT1 (|2010D was a breakthrough for the 
physics of neutron stars. If a more massive pulsar was de- 
tected in the future, the constraints would become even 
more severe. However, pulsar mass measurements give 
always only a lower limit for the actual maximum mass. 
Contrary, with the black hole formation one has direct 
access to the (proto-neutron star) maximum mass. Fur- 
thermore, if we would know the real maximum mass of 
cold neutron stars, then one could use the neutrino signal 
to constrain the temperature dependence of the EOS. On 
the other hand, the correlation of Fig. [17] has also a lot 
of predictive power if it is combined with the mass con- 
straint of pulsar PSR J1614-2230. All EOS with black 
crosses to the left of this line are not compatible with 
the observational data of PSR J1614-2230, because their 
maximum masses of cold neutron stars are too low. Thus 
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Fig. 17. — Time fo r the onset of collapse to a b lack hole for the 
40 Mq progenitor of Wooslcy and Weaver (1995) with respect to 
different configurations of maximum gravitational masses for the 
investigated EOS. "cold": static cold neutron stars at T = 0.1 MeV 
and beta-equilibrium without neutrinos, "s = 4": static proto- 
ncutron stars at a constant entropy per baryon of s = 4 kg and 
beta-equilibrium without neutrinos, "sim" : the maximum gravita- 
tional masses found in the simulations. The vertic al line shows the 
mass measurement of pulsar PSR J1614-2230 by Dcmorcst ct al. 

ipmToT i. 

we expect that the neutrino signal of such a black hole 
formation event from a massive 40 Mq progenitor will 
not be much shorter than 0.5 s. On the other hand, if 
the signal would end at t p h <C 0.5 s this would indicate 
the missing of important physics in the current models. 

Finally, we have to make a remark of caution for the 
results of this section: for such high temperatures as 
encountered here, e xtensions of the EOS m ay become 
important, see, e.g.. lOertel fc Fantinal (|2010l ). Tempera- 
tures of 100 MeV are on the order of the pion mass, so 
a pion gas would form, and also muons and baryon reso- 
nances could give an important contribution to the EOS. 
In addition one has to be cautious using a non-relativistic 
EOS at such extreme conditions. It is not warranted that 
the EOS remains causal, i.e., that the speed of sound 
stays below the speed of light. We find that LS becomes 
acausal during the collapse to a black hole. However, we 
do not expect that our results will change qualitatively 
by these aspects, but still they should be studied more 
elaborately with extended EOS in the future. 

5. SUMMARY AND CONCLUSIONS 

We apply the new RMF EOS of t he statistical model of 
iHempel fe Schaffner-Biehchl (|2010D (HS) to core-collapse 
supernova simulations of 15 and 40 M Q progenitors. Our 
core-collapse model is based on general relativistic ra- 
diation hydrodynamics that employs three-flavor Boltz- 
mann neutrino transport. The HS EOS include up to 
8000 nuclei , base d on experimentally known masses of 
lAudi et all (|2003l ) and different theoretical nuclear struc- 
ture calculations. Three different parameterizations of 
RMF interactions are applied, TM1, TMA, and FSUgold 
which govern the high-density behavior of the EOS. Re- 
garding the nuclear matter properties and mass-radius 
relation of cold neutron stars, none of the HS EOS used 
here was able to fulfill all of the quoted experimental and 
observational constraints. The same h o lds for the com- 
monly used E OS from H. IShen et all dl998ai) (STOS) 
and the EOS of lLattimer fe Swestvl ( |1991ft (LS) with an 
incompressibility of 180 MeV (LS (180)). Furthermore, 



we also include the LS EOS with an incompressibility of 
220 MeV (LS (220)), which is in reasonable agreement 
with the considered constraints. The HS EOS tables are 
available for use and can be accessed from the home page 
given in footnote [T] Also routines are available which 
allow to calculate the abundances of all considered nu- 
clei. Additional EOS tables for the NL3 (|Lalazissis et all 
[19971 ) and DD2 (|Typel et al.|[2010l ) RMF parameteriza- 
tions can also be downloaded, which were not included 
in the present study. We remark that DD2 is in good 
agreement with the quoted experimental and observa- 
tional constraints. In the future, the EOS tables will 
also be provided on the CompOSE(db) online platform 
of the CompStar research networking program. 

We identify the most important differences between 
the HS EOS with the TM1 parameterization, the STOS 
EOS which is also based on TM1, and the LS (180) EOS 
by simulating the evolution of the 15 M progenitor. 
For the first time, two EOS are available which have 
identical nucleon interactions but a different description 
of the low-density phase where nuclei exist. In general 
one of the big advantages of the HS EOS, in compari- 
son to LS and STOS, is the detailed information of the 
nuclear composition, the precise description of nuclei in- 
cluding shell effects, and the consideration of additional 
light nuclei. Regarding the first two aspects, the low- 
temperature HS EOS leads to a smooth transition from 
NSE to non-NSE at temperatures of about 0.5 MeV. 
There, HS matches the ideal gas of 56 Fe or 56 Ni (de- 
pending on the proton-to-baryon ratio) by construction, 
while LS and STOS produce artificial jumps, e.g., in in- 
ternal energy. Also the agreement with the composition 
and the thermodynamic state of the progenitor was most 
satisfactory with HS. In the same sense, the HS EOS is 
well suited to be connected to explosive nucleosynthesis 
calculations, as it naturally gives a realistic seed distri- 
bution. Obviously, this is true for all NSE models which 
employ experimental nuclear binding energies, but the 
HS EOS tables also cover all densities and temperatures 
relevant for supernova simulations, which is not typical 
for most of the NSE models. 

During the collapse we observe that HS (TM1) gives 
significantly smaller nuclei than STOS (TM1), leading to 
a stronger deleptonization and a lower electron fraction 
in the core. Contrary, the lower free proton fraction in 
the outer layers using HS caused less electron captures 
there. This illustrates impressively the sensitivity of the 
core-collapse simulation to slight changes of the nuclear 
composition. The electron captures in HS are less en- 
ergetic, partly due to the use of real nucleon masses in 
the EOS, which therefore contain the proton-to-neutron 
rest-mass difference. During the collapse HS leads there- 
fore to less entropy generation in the center compared 
with STOS. 

The collapse phase is dominated by heavy nuclei. We 
investigated the nuclear distributions and identified a vis- 
ible effect of neutron magic shells in the nuclear compo- 
sition, leading, e.g., to several peaks in the mass distri- 
butions, which cannot be captured by the SNA. Inter- 
estingly, closed neutron shells also have an influence on 
the electron fraction profiles. In the innermost zones the 
nuclear distributions become very broad during collapse, 
due to the large compression heating. It would be a very 
appealing extension of the HS EOS model to calculate 
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weak reaction rates with the nuclear distributions of the 
EOS. Here, we only considered the avera ge heavy nucleu s 
in combination with the reaction rates of lBruennl (1985). 

The appearance of additional light nuclei, which so far 
were not implemented in EOS tables for core-collapse 
supernovae, was investigated in detail. The formation 
of weakly bound light nuclei is driven by entropy. They 
are formed in the supernova by an interplay of density 
and temperature, where the temperature has to be high 
enough to suppress the role of differences in binding en- 
ergies. We find that during the early collapse phase, 
the commonly used approximation to consider only al- 
pha particles of all possible light nuclei is acceptable. 
However, already a few milliseconds before bounce tem- 
peratures are sufficiently high so that additional light 
nuclei can be even more abundant than alpha particles. 
When the shock runs through the infalling matter the 
well-known dissociation of heavy nuclei occurs. In the 
present article we further show that a significant frac- 
tion of light nuclei, e.g., 2 H, 3 H, 3 He, and Li, forms in 
the hot dissociated matter between the shock and the 
core of the proto-neutron star. These light nuclei remain 
there also in the post-bounce phase. For the inner lay- 
ers of the proto-neutron star envelope, where matter is 
very asymmetric, tritium gives the largest contribution. 
In the outer layers below the shock, deuterium plays 
the most important role. We remark that these results 
for light nuclei obtained with the HS EOS are in qual- 
itative agreement with e xpectat i ons fr om t he quantum 
many- bod y approaches o f Ropke (20fl|) and lTvpel et al.1 
(201 0h fseelHempel et all (120111) 1 and the previous study 
of lSumivoshi fc Ropkd (|2008D . 

In the present article, we simplified the neutrino reac- 
tions with light nuclei by treating all of them as alpha 
particles, i.e., we only considered elastic neutral-current 
reactions. However, light nuclei have low binding ener- 
gies, except alpha particles. Because of the lower binding 
energies we expect that weak reactions with the other 
light nuclei will be different compared to the weak reac- 
tions with alpha particles. Some aspects may resemble 
the reactions with free nucleons. For example, also elec- 
tron captures or break-up reactions of the deuteron are 
likely to occur. Until 50 ms post bounce the abundances 
of deuterons and tritons close to the shock are compara- 
ble to the ones of free protons. Thus, there may be an 
influence of light nuclei on neutrino heating during the 
early post-bounce phase. However, the fractions of light 
nuclei in the heating region decrease in the subsequent 
evolution. The largest fractions of light nuclei are then 
found below the neutrinospheres. Thus for the later post- 
bounce evolution cooling may be affected rather than 
heating. More detailed studies are required to answer to 
which extent light nuclei can contribute to the neutrino 
heating and cooling and possibly on the explosion dy- 
namics. We see the consistent implementation of weak 
processes with light nuclei in the neutrino transport of 
core-collapse supernova simulations as an interesting task 
for future studies. 

Regarding the hydrodynamic behavior of the new 
HS (TM1) EOS, we observed that it leads to a more 
compact proto-neutron star with slightly higher temper- 
atures than with STOS (TM1). This is a result of the 
different description of the non-uniform nuclear matter 
phase, giving the appearance of light nuclei and smaller 



heavy nuclei, which in turn leads to lower central electron 
fractions in the simulations. During the post-bounce evo- 
lution, the differences between HS (TM1), LS (180), and 
STOS (TM1) reduce. However, the much softer LS (180) 
EOS has significantly higher central temperatures and 
densities, and hence the standing accretion shock con- 
tracts faster. We arrive at the general conclusion that 
the model for low-density nuclear matter, where nuclei 
exist, and the implemented nuclear degrees of freedom 
are as important as the description of the nuclear inter- 
actions around and above saturation density. 

We also compare the evolution of neutrino luminosi- 
ties and mean energies for the different EOS, for the 15 
M Q model. The soft LS (180) EOS leads to high neu- 
trino fluxes and mean energies. Although STOS (TM1) 
and HS (TM1) are based on the same RMF parame- 
terization, the larger deleptonization for HS results in 
larger neutrino fluxes and mean energies than STOS, 
emphasizing our statement above. The neutrino signa- 
tures are related to the evolution of the neutrinospheres 
which are located at the proto-neutron star surface. For 
LS (180) the proto-neutron star contracts fastest with 
highest temperatures and Y e obtained, while STOS and 
HS have significantly longer contraction times especially 
after 300 ms post bounce. However, the differences ob- 
tained in the neutrino luminosities and mean energies are 
rather small. Therefore for such intermediate-mass pro- 
genitors, very precise neutrino observations of a possible 
future galactic supernova were required to constrain the 
nuclear EOS further. Note, this statement remains only 
valid if EOS effects are not affected significantly by the 
explosion, which is expected to occur for such a progen- 
itor. 

Pronounced EOS effects are found for the neutrino 
signal of the core-collapse supernovae of a 40 M Q pro- 
genitor. We confirmed the results of other studies 
that /i/r- neutrinos carry an important signature of the 
EOS, namely, the temperatures reached inside the proto- 
neutron star, which are connected to the stiffness around 
saturation density. However, the most prominent EOS 
aspect related to the neutrino signal of such massive pro- 
genitors is the time until black hole formation. We found 
that it is dictated by the global behavior of the EOS 
and not by single EOS parameters. We showed that the 
correlation between time until black hole formation and 
maximum mass of cold neutron stars, found in previ- 
ous EOS studies, does not hold for the extended set of 
EOS investigated here. Using the two non-relativistic LS 
EOS, the collapse to the black hole takes place system- 
atically earlier than with the RMF EOS, even if the LS 
EOS have larger maximum masses for cold neutron stars. 
This interesting result was explained by the different re- 
sponse to finite temperature. For the RMF EOS a signif- 
icantly larger increase of the maximum masses at finite 
entropy was observed than for the LS EOS. The states at 
the beginning of black hole formation in the simulations 
could be approximated by hydro-static configurations of 
proto-neutron stars at a constant entropy per baryon of 
s = 4 kfl. The resulting maximum masses agreed with 
the results from the simulations. We found a new mono- 
tonic correlation between the time until black hole for- 
mation and the maximum masses of these proto-neutron 
star configurations. From th is correlation and the max- 
imum mass constraint from iDemorest et al.l (|2010f) , we 
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predict that the neutrino signal of a core-collapse su- 
pernova of a 40 M progenitor should last at least for 
0.5 s starting at bounce, until the black hole forms. If 
such an event was measured in the future, the correlation 
would provide the maximum mass of proto-neutron stars 
at constant entropy. This would represent a significant 
constraint for the nuclear EOS at finite entropy, comple- 
mentary to observations of cold neutron stars. In combi- 
nation with the existing knowledge about the maximum 
mass of cold neutron stars this would give information 
about the temperature dependency of the EOS. 

The EOS can also affect the explosion dynamics. In 
general, soft EOS lead to a stronger initial shock acceler- 
ation where higher matter velocities are o btained in com - 
parison to stiff EOS (jBaron et al.lll985allb1 : lB~ruennlll989[ ). 
Therefore soft EOS result in a more optimistic situation 
with respect to the initial energetics at shock breakout 
shortly after bounce. This was confirmed in the present 
study, where the highest matter velocities for the bounce 
shock were encountered for the extremely soft LS EOS, 
and the lowest for the stiff TM1 EOS. Furthermore, al- 
though based on the same RMF parameterization TM1, 
STOS and HS have different shock energetics at shock 
breakout shortly after core bounce, due to the different 
description of nuclear matter below saturation density, 
and the resulting different deleptonization during the col- 
lapse. However, the situation becomes more similar after 
the bounce shocks have stalled. The initial dynamical 
advantage of a soft EOS is lost with time. It neverthe- 
less leads to a more extended heating region between the 
neutrinospheres and the standing bounce shock. The in- 



fl uence of the EOS on ex plosions has been investigated 
in iMarek &: Jankal (|2009l ) based on simulations in axial 
symmetry. It would be an interesting follow-up project 
to perform a similar EOS comparison, as was done here, 
in multi-dimensional simulations, to further tackle the 
explosion mechanism of core-collapse supernovae. 
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